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Abstract 

This paper is extended and broadly generalized version of earlier pub- 
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on some problems which were left unsolved or just mentioned in Phyiscs Re- 
ports 298, 251 (1998). Applications of the obtained results include (but 
not limited to)polymer physics, classical and quantum chaos, fractional Hall 
effect, dynamics of textures in liquid crystals and 2+1 Einsteinian gravity. 
The paper presents some self contained excerpts from mathematics (not dis- 
cussed so far in physics literature) to facilitate the uninterrupted reading of 
the manuscript. 

PACS number(s):02.50.-r,05.40.+j,61.41.+e 



*375 H.L. Hunter Laboratories ,Clemson University, Clemson ,SC 29634-1905 , USA 



1 



I. Introduction 



In our recent short paper [1] we had provided new derivation of Spitzer's 
law [2] for Brownian windings . This problem is of some importance in poly- 
mer physics since it allows to define the concept of entanglement in precise 
mathematical terms. It is also of relewance to other disciplines , e.g. 2+1 
gravity , quantum Hall effect, etc. Due to format limitations of our paperfl] 
many important topics were either briefly mentioned or not discussed at all. 
This paper serves to correct this deficiency. If the reader will have patience 
to go over the entire manuscript , he or she will recognize, that, still, this 
much longer paper is just an introduction into the very rich field of various 
topological problems which are waiting for development. In particular, in 
this work we had only briefly mentioned the hyperbolic groups (in section 
6) while this has become a very large area of research in mathematics [3] 
.Study of random walks on hyperbolic groups and their associated graphs 
reveals that, for example, the notion of ultrametricity which is used in the 
theory of spin glasses[4] is just a special case of hyperbolicity.Hence, the 
static and dynamic problems associated with spin glasses, glasses in general, 
and all other cases which involve replicas, could be studied without replicas 
by properly formulating the ultrametric problems within the scope of the 
hyperbolic group methods. The steps in this direction were already made by 
mathematicians [5] but large amount of work is still waiting for development. 
Also, the issues related to symbolic dynamic and finite state automata inter- 
pretation of the results presented in this paper were left completely outside 
the scope of this paper. Interested reader can ,in part, correct this deficiency 
by reading recent review article by Adler[6]. 

This paper had been also motivated to some extent by the recent papers 
in physics literature where one and the same mathematical object is being 
treated by completely different methods. Specifically, in Ref.[7] the authors 
study the extension of the results of Aharonov-Bohm [8] to the case of two 
magnetic fluxes while in Ref.[9] quantum mechanics in the rectangular bil- 
liard with single pointlike scatterer was studied. Subsequently, the results 
of the last paper were improved [10] and generalized to many scatterers in 
Ref.fll]. To this list of papers ,one may add Ref. [12], etc., while in Ref.[13] the 
random walks on the disordered cellular networks were studied. In this pa- 
per it is shown, that all these problems are interrelated and ,therefore, could 
be treated using the same mathematical methods. Moreover, other problems 
could be studied by these methods as well and some of these problems are 
discussed in section 7 . In particular, by analogy with the Aharonov-Bohm 
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problem in two dimensions which had become a sort of Hydrogen atom-like 
problem to all problems which involve physical processes in the multiply con- 
nected spaces, in three dimensions one can think of a problem about random 
walks in the presence of nontrivial knots mentioned already in our earlier 
work, Ref[14]. Such problem may have at least some biological significance 
since it is important to know to what extent diffusion processes are being 
affected by the topology of ,say, knotted DNA. In this work, we use the 
same methods as that developed for the treatment of random walks on fig- 
ure eight in order to provide some solutions to the problem of random walks 
in the presence of knots. Alternative methods of solution of this problem 
are discussed in Ref.[15]. We also reconsider the dynamics of area-preserving 
toral homeo(diffeo)morphisms[16] from the point of view of random walks 
on the Teichmiiller modular group. Such reconsideration is useful for several 
reasons. First, it allows to discuss new ways of classifying chaotic and non- 
chaotic systems through the associated with them 3-manifolds. Second, it is 
very useful for the description of dynamics of textures in liquid crystals and 
particles in 2+1 gravity . In this paper,naturally, above topics are studied 
only in passing so that interested reader is encouraged to read our recent 
papers [17] for better understanding of all aspects of these problems. In 
this work we demonstrate how the figure 8 problem is related to the so- 
lution of David Hilbert's 21st problem( in section 5). In the same section 
we discuss how this problem is connected with the theory of exactly inte- 
grable systems, Knizhnik-Zamolodchikov equations, quantum groups , theory 
of elasticity. We would like to emphasize that our presentation is not just 
a review of literature: each section contains worked out examples which are 
original in nature. At the same time, some mathematical facts are presented 
along with these developments to facilitate the uninterrupted reading of the 
manuscript. 

In section 2 some old results related to random walks on once punctured 
plane are revisited from the point of view of the results of Ito and McK- 
ean[18], Lyons and McKean[19] and McKean and Sullivan[20]. The rest of 
the paper could be considered as various attempts to extend ,to reanalyze 
and to generalize their results to the case of multiple punctures .The guiding 
principle in doing this is associated with various aspects of random walks 
on groups , mainly free groups. The comprehensive mathematically rigorous 
and physically readable review of this subject can be found in Ref.[21]. 
In our work only very small portion of the available mathematical results 
is utilized. In section 3 we introduce some concepts related to groups and 
group presentations in order to use them in section 4 devoted to finding 
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the function which connects the multiply connected trice punctured sphere 
- the space in which the random walk on figure 8 is taking place-to the 
simply connected universal covering Riemann surface. Although in case of 
trice punctured sphere such function can be found generalization of the 
methods by which it was obtained to more than three punctures represents 
a formidable task associated with mathematical methods which so far had 
been used only in string theory. We discuss these methods in the same section 
and later,in section 7 , provide an introduction to the completely different 
approach to the whole problem based on mathematical ideas of Grotendieck 
and Belyi [22]. The results of section 4 are being put in a broader context of 
Riemann-Hilbert problem and Knizhnik-Zamolodchikov equations in section 
5. This section also clarifies the connections with the gauge field-theoretic 
treatments of the same problems popular in physics literature [14]. In sec- 
tion 6 we develop some results of Kasteleyn[23] related to random walks on 
free groups. Based on this development, we reobtain the results of Lyons and 
McKean for transience/recurrence of random walks on the trice punctured 
sphere. We also provide in this section a definition of hyperbolic groups and 
hyperbolic spaces and show that the ultrametric spaces (used in the theory 
of spin glasses[4]) are just special cases of hyperbolic spaces. The results 
obtained in this section are further developed in section 7. In this section, 
solely devoted to applications,we discuss several topics. In particular, we dis- 
cuss random walks in the presence of knots. We demonstrate that the tran- 
sience/recurrence of random walks in the presence of a knot depends upon its 
nature. For example, random walk on figure 8 is related to the random walk 
in the complement of the figure 8 knot and , hence, is transient while random 
walk in the complement of the trefoil knot(the simplest among torus-type 
knots) is recurrent. Next, we discuss the chaotic v.s. regular behavior of 
such walks by analyzing their behavior with help of number-theoretic meth- 
ods, by using the theory of Teichmuller spaces and by employing some facts 
from the theory of 3-manifolds. These methods provide new efficient ways 
of description of chaotic systems and are potentially useful for treatments 
of dynamics of liquid crystals and 2+1 gravity. Naturally, these topics are 
treated only superficially in this paper, but ,we hope, that many illustrations 
which we supply along the way, should be helpful for better understanding 
and appreciation of much deeper results available in mathematics [24] and 
used in our earlier works[17] for treatment of these physical problems. This 
section is finished with several additional applications. For instance, use 
of Markov triples , known in the number theory for more than a century 
[25] along with the results of Grotendieck and Belyi [22] provides us with an 
opportunity to introduce an alternative explanation of the fractional Hall 
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effect. These results are also complementary to that presented in section 
4. We also discuss some connections between the results the theory of free 
groups and the theory of braids . Such discussion reveals advantages and 
disadvantages of either methods for studying evolution of random walks in 
multiply connected spaces. 

II. Some results about random walks on once punctured 
plane 

Study of random walks on once punctured plane is closely related to 
the Aharonov-Bohm (A-B) effect in quantum mechanics [ 8]. Although this 
problem was studied very thoroughly, here we would like to emphasize some 
of its aspects which, to our knowledge, had received much less attention. 
In particular,we would like to investigate how the description of random 
walks on once punctured plane R 2 — is related to the description of the 
random walks on the circle S 1 .The last problem was studied in great detail by 
Schulman [26 ,27] and is considered to be the benchmark problem associated 
with quantization in multiply connected spaces [28 ]. 

Although it is intuitively clear that both problems are connected , to 
establish such connection mathematically is nontrivial.Let us begin with 
the observation that in both cases we are dealing with multiply connected 
spaces. The theory of differential equations on the punctured plane was ini- 
tiated already by Poincare and is known in the literature as the theory of 
Fuchsian equations [29,30 ]. According to this theory ,it is essential to find 
a simply connected (covering )space which is associated with multiply con- 
nected (base) space. Then, the related problem is solved by lifting it to this 
covering (simply connected )space and then,the solution is projected down 
to the base space. To find a transformation from the base space space M 
(multiply connected) to the covering space M (simply connected )could be 
already a very difficult (and even unsolvable !) problem .Accordingly,some 
simplifications may be required in order to obtain at least a partial informa- 
tion. The circle S 1 and the punctured plane R 2 — represent a good examples 
of use of general methods outlined above and are ideally suited for study 
of simplifications. Fortunately ,both problems could be solved completely 
so that they could be considered as having the same status in the theory 
of diffusion/quantum mechanical problems in multiply connected spaces as 
the Hydrogen atom problem in standard quantum mechanics. 

The universal covering space M for the circle problem is real line R so 
that S 1 = R/Z where Z is the set of all integers .For the random walk of 
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N effective steps on R we expect our distribution (Green's) function to be 
transitionally invariant , i.e. 

G%(x',x) = G%(x'-x), (2.1) 

where is the Green's function (end-to-end distribution function) for the 
random walk on R and x(x') is the initial (final )position of the walk in this 
space. Surely, Eq.(2.1) can be equivalently rewritten as 

G^(x' + 2irn,x + 2Trn) = G%(x',x), (2.2) 

where n = 0,±l,±2,... The above example can be easily generalized now. 

Let 7 be some generator of motion(s) on M ,7 £ T, where T is some 
group of motions on M , then,in general, we should require 

G§( 7 x > , 1 x) = G§(x',x). (2.3) 

since this result holds for Eq.s(2.1) and(2.2). The above equation can be 
equivalently rewritten as 

G§( 1 x',x) = G§(x', 1 -^)- (2-4) 

This follows easily if we replace x by 7~ 1 x in Eq.(2.3). In the case of 
quantum mechanics the transition from M to M is made via 

Gft(x',x) = ^2p( 1 )G§( 1 x',x) (2.5) 
where p(x) is the phase factor , 

p( 7 ) = e l5(7) (2.6) 
which is subject to the cocycle constraint [31 ,32 ] 

*(77) = %) + %'). (2.7) 

This constraint follows easily from the composition law (Markov property) 
for propagators. 

In the case of punctured plane the phase factor 5 is actually responsible 
for changes in statistics and/or for particle- hole interactions. In the case of 
S 1 the phase S = 0(or 5 = it) is associated with bosons (or fermions) [26]. There 
is yet another interpretation of the phase factor. It comes from the polymer 
physics .In Ref.[14] we had demonstrated that the description of the random 
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walks on S 1 can be associated with the description of random walks (poly- 
mers)which are confined between the parallel plates (or planes)with which 
these walks are allowed to interact. The interactions are effectively introduc- 
ing statistics so that , depending on their strength,the polymer chain acts 
as bosonic or fermionic quantum particle or as a particle with fractional 
statistics. 

The result given by Eq.(2.5) still deserves some additional comments. Let 
us discuss the case of a simple random walk so that we can put £(7) = 
0. Already in Eq.(2.2) we had introduced the fundamental domain D: 0< 
x < 2tt, so that the translations producing R are created by the successive 
use of 7 applied to some xG D.This simple result should hold in more general 
cases too. That is,we expect that the manifold M should be expressible as 
a quotient M = M/T .More intuitively,M should be related tosome funda- 
mental domain in M so that the whole M is being covered (without gaps!) 
by translations of the fundamental domain with help of the group elements 
7 G r. Consider , therefore, instead of Eq.(2.5),its modification 

G§(x',x) = ^2G§( 1 x',x). (2.8) 

Let 7W G T. We would like to demonstrate now that 

G§ d (1) x', 7 «x) = G§ Or', x) = G§ (x' - x). (2.9) 
Indeed, using Eq.(2.8),we obtain 

G§ ( 7 «x',7 (1) x) = ( 77 Wf', 7 Wf). (2.10) 

Using Eq.(2.4) ,we obtain as well 

E °N ( (7 (1) ) _1 77 {1) i', x) = Y,Gn (7^, 5) (2-11) 

since (7^^) 1 77W G 7 G T.This proves Eq.(2.9). Using this result ,it is 
possible to prove the following 

Theorem. 2.1. Let G^(x',x) be the distribution function for the random 
walk on the covering space M which possesses property reflected in Eq.(2.9) 
,then the distribution function on the base space M can be obtained from 
G^J through simple identification: 

G§(x',x) = G^(x',x) (2.12) 



7 



where x' and x are inverse images of x and x' under the natural projection: 
M -> M = M/T. 

Proof. Please, consult Rcf. [28] . 

Remark. 2. 2. Using Eq.(2.9),it is clear , that it is sufficient that x and 
x 1 G D.This means, that if we are able to find solution for the random walk 
problem in the covering space M ,we can obtain solution of the random walk 
problem in the base space M without additional complications. 

Let us illustrate the above general statements on example of the random 
walk on R 2 — 0. Using dimensionless units we obtain the following diffusion 
equation , 

°l-l(*_ + *\ f (213) 

8t "4 V0s 2 03/V { ' 

where the factor | (the diffusion coefficient) is chosen for further convenience 
and / = f(x, y; t) with t ^ N. It is convenient to rewrite Eq.(2.13) in terms 
of complex variables:^ = x + iy,z = x — iy. Simple calculation produces then 

df d 2 

■s = ear,f- (2 - 14) 

Since this equation is not defined for z = z = 0, we would like to introduce 
new complex variable w via z = exp{w}. Unlike z,the complex variable 
w = u + iv is defined on the entire complex w— plane. In terms of u;— variable 
Eq.(2.14) can be rewritten as 

V = e ~2u^ J 

dt dwdw 

Obtained equation describes Brownian motion on a simply connected cov- 
ering space M which is the Riemann surface of the logarithmic function, i.e. 

w = \n\z\ +i(argz + 2nn),n = 0, ±1,±2, (2.16) 

and is known to be simply connected. The result just obtained coincides with 
that discussed in the book by Ito and McKean,e.g.see page 280 of Ref.[18 
].Ito and McKean argue (without demonstration) that the diffusion Eq.(2.15) 
can be converted into the form given by Eq.(2.13) by replacing Brownian 
time t by another , actually, random time T. 
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The idea of the proof goes back to Paul Levi [33] and was refined by 
others [34]. Let us discuss in some detail how this can actually be done. To 
this purpose , let us rewrite Eq.(2.15) in the following equivalent form 

^> -I (2 17) 

dt e ~4\du 2 0wV [ ' 

Now let us select new time T by requiring 

^ = eM-MT)} 

or,equivalently, 

T 

t = j drexp{2u(r)}. (2.18) 
o 

This then will bring Eq.(2.17) back to the form of Eq.(2.13) as anticipated. 
Eq.(2.18) coincides with that presented on page 280 of Ito and McKean [18 
] where it was given without derivation. In spite of its simple form , thus 
transformed diffusion equation will depend upon the random time T since 
u(t) is Brownian motion. Accordinglythe observables calculated with help of 
such diffusion equation will require an extra averaging .This may sometimes 
be inconvenient in practical calculations .To recognize this difficulty we shall 
discuss a generic example which exhibits all the features involved. Before 
doing so , several remarks are in order. 

First, it is clear ,that in more general case of multiple punctures one 
has to find an analogue of Eq.(2.16) ,i.e. to find some function z = &(w) 
such that w— plane (actually Riemann surface) is going to be an universal 
covering space M which is simply connected by construction. 

Second, to find such function is not an easy task and may not be possible 
in general, e.g. sections 4 and 7 below. Nevertheless, if such function can be 
found,then, by analogy with Eq.(2.15) ,we will be able to write 

8/ ^--ViU^/ 



$ (w] 



\ du 2 dv 2 J 



dt 4 

and, accordingly, the new time can be defined as 

T 



t= dr 



2 



$ (w(t)) (2.20) 



o 
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so that the diffusion equation ,Eq.(2.19) , again acquires its canonical form 
given by Eq.(2.13) in accord with general result by Paul Levi [33]. 

Let us now return back to our example of once punctured plane. In 
this case, the explicit form of the end-to- end distribution function is well 
known, e.g. see Refs. [14, 18,35 ] ,and in polar system of coordinates it is given 
by 

G(ri,r 2 ,A0;t) = — exp{-^±^} £ e* mAe I m (z) (2.21) 

m=— oo 

, where A9 = 02—Oi,z = 2r\T2jt and I m (z) = I- m (z) is the modified Bessel's 
function. The above distribution function can be used either for study of the 
radial or the angular distributions or both . Suppose,we are interested in 
the angular distribution function only .Then, using Eq.(2.21),it is convenient 
to introduce the normalized distribution function defined according to the 
following prescription : 

^-W^-mjL*" 1 ""- (2 ' 22) 

The Fourier transform of such defined distribution function can now be 
obtained in a standard way as 

oo 

F(z,a) = [ dA9 e~ iaAe F(z, AO) = . (2.23) 

J hyz) 

— oo 

Let us now choose r 2 = rs/i + ri.This choice is motivated by known scaling 
properties of Brownian motion [34]. For large t one obtains :z ~ 2r\f j \ft. 
For fixed f and r\ and t — > oo one surely expects z — > 0. Using known 
asymptotic expansion of I\ a \(z) for small z's in Eq.(2.23) allows us to obtain 
closed form analytic expression for F(z,a) in this limit as 

led 

F(z,a) ~ exp{--^-lnt}. (2.24) 

The inverse Fourier transform of Eq.(2.24 ) produces famous Cauchy-type 
distribution for AO 

1A6 1 1 

P(x = )dx = ^dx (2.25) 

v lnt ; vrl + x 2 y ' 

obtained originally by Spitzer in 1958 [2] in much more complicated way. 
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Let us now reproduce this result with help of the random time change 
discussed above. To this purpose, let us begin with the classical Lagrangian 
L for the fictitious Brownian particle traveling in the punctured plane. By 
introducing polar system of coordinates our L can be written as 

L[r] = l(r 2 + rV). (2.26) 

The path integral for such particle can be written now in a usual way as 

t 

G( ri ,r 2 ,A9;t) = J D[r[r],9[T]] exp{- J cItL[t}} (2.27) 

o 

where the limits of path integration are not written explicitly. Following the 
same line of arguments as that which had lead us from Eq.(2.13) to (2.15) 
let us represent r— variable through x— variable defined by r = expx. Surely, 
for — oo < x < oo,we have < r < oo as required. The Lagrangian L is now 
being converted to 

L = ^ e 2x (x 2 + 9 2 ) (2.28) 

which, in view of Eq.(2.27) , should have its analogue in Eq.(2. 15). Consider 
now the change of time variable inside of the path integral in Eq. (2. 27). That 
is we have r = r(T) so that ,for example, 

dx dT , , , , / dx \ 2 , m dT ( dx \ 2 

— — and,accordmgly, dr — = dT — — . (2.29) 



dT dr ' \dr J dr \dT 

This then requires us to make a choice for ^ which is identical with that 
made in Eq.(2.18) so that in terms of new time T our path integral,Eq.(2.27), 
represents, indeed ,a free diffusion process in which x and 9 are two inde- 
pendent Brownian motions in accord with our previous calculations. Let us 
now calculate the generating function F(z, a) which was defined in Eq.(2.23) 
but ,this time,with help of the propagator Eq.(2.27) (with time rescaled).We 
have 

I iajdr6{T)\ 

F(z,a) = (e o \ (2.30) 

where (...) denotes functional integral averaging with help of the properly 
normalized propagator , Eq. (2. 27). The Gaussian integration can be trivially 
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performed now with the result : 

oo x 2 (T) T 



F(z,a)= J dT J D[x(T)]5(t- J dr exp{2x(r)}) exp{-^ J dr(i 2 + a 2 )}. 



xi (o) 



(2.31) 



In arriving at this result we took into account the constraint ,Eq.(2.18),so 
that the final answer is written in terms of the unrescaled time. The con- 
straints of this sort were discussed before in Refs [36, 37]. The presence of 
constraint complicates calculations somehow. To by- pass this difficulty,the 
Laplace transform of Eq.(2.31) can be taken with the result 

oo x(T) T 

F s (z,a) = J dT j £>[x(r)]exp{-i J dr(x 2 + a 2 + 2s exp{2x(r)})} 



s(0) 



(2.32) 



where the Laplace variable s is conjugate to time r .The corresponding 
Schrodinger-like equation can be now easily written as 



sexp{2x} - -a 2 ) 



F s (z, a) = 5(x — x ). 



(2.33) 



By using again r-variable(instead of x) we obtain instead of Eq.(2.33) the 
following equivalent equation (x / x') 



sr 



d 2 



+ r- 



dr 2 dr 



a 2 ) 



F s (z,a)=0 



(2.34) 



in which one can easily recognize the equation for the modified Bessel's 
function(s) I a and K a . The Green's function can be constructed with help 
of these functions in a standard way,e.g.see Ref.[18] ,with the result : 



F s (z,a) 



7| a |(v / 2sr 1 )K| a |( v / 2sr 2 ) 
K\ a \{VTsn)I H {VTsr 2 ) 



if r\ < r 2 

if T2 < T\ 



(2.35) 



With help of this function the inverse Laplace transform now can be per- 
formed with help of Eq.(56) of Cr.5 of Bateman and Erdelyi [38] resulting 
in 



2 2 

F(z,a) = -exp{- r -l±^}I H (z) 



(2.36) 
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It is clear ,that this result coincides with earlier obtained Eq.(2.23) provided 
that it is properly normalized. 

Prom the above derivation it is clear that the entire distribution func- 
tion, Eq. (2. 21), can also be obtained. By repeating the same steps as we have 
used to arrive at Eq.(2.33) ,we obtain now 



rnl i , d 2 d\ 



G(x,x ,9) = S(x-x)S(9-9). (2.37) 



This equation is defined in the fundamental domain in the covering space M 
which is just strip of finite width 2ir in 9 direction and infinite -oo < x < oo in 
x direction. Naturallythis equation coincides with Eq.(2.17) (if the Laplace 
transform of Eq.(2.17) is performed). It is remarkable that the radial and the 
angular Brownian motions are completely decoupled in terms of variables 
which we are using. 

If,as before ,we introduce r-variable via r=expx ,then Eq.(2.37) becomes 
again Bessel-type equation for the corresponding Green's function. The fun- 
damental domain D under such change will be converted from an infinite 
strip to the whole r — 9 plane with cut along the positive semiaxis.In this case 
we observe that Eq.(2.21) is in accord with Eq.(2.12) as required. Evidentlywe 
could as well use the strip of width tt ,then the strip is converted into the 
upper half plane. The Riemann surface of the logarithm can be made by 
properly gluing together either the stuck of cut r — 9 planes or half planes 
.That is the Riemann surface can be described (or constructed) in several 
ways [38]. The half plane description is especially useful since it allows to 
obtain the Cauchy-type distribution at once [34]. At the same time,the de- 
scription in terms of r — 9 variables leads to the undesirable coupling between 
the radial and the angular Brownian motions. Indeed,the Green's function 
for the circle S' 1 of radius r is known [26,31] to be expressed as 



1 2 

G(A9,t) = — V e imA V^ . (2.38) 
2tt z — ' 



For r\ = r2 and z>l use of the asymptotic estimate [40] 

J »=vfe (1 + i +0 ( 2 " 2 » ex ^-£ +0(2 " 2 » 

allows us to bring Eq.(2.21) into the form 



G( n = r 2 = r, A9; t) ~ \/ ^ e imM e~^ . (2.39) 

m=— oo 
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Hence,in terms of usual polar coordinates, it is impossible to disentangle the 
angular and the radial parts of Brownian motion for arbitrary values of r 
and t . At the same time, depending upon the questions being asked, it may 
be sufficient to use ,say,the angular part, e.g. Eq.(2.38) only .Moreover,the 
obtained results can be generalized to d-dimensions.In this case the Green's 
function analogue of 5 1 problem is known to be [18] 

oo 

G(e 1 e 2 ;t) = J2^ tJn E s n(0i)^(0 2 ) (2.40) 

n=0 l<m(n) 

with 7„ = \n{n + d — 2) and S l n (9) obeying the generalized spherical har- 
monics equation 

l -^S l n = ln S l n . (2.41) 

The above Green's function should be compared with that for R 2 — .As 
results of Ref.[18] indicate,the reduction analogous to that given by Eq.(2.39) 
is not possible for d>2. 

In the case of random walks the most commonly asked questions are 
associated with the probability of returning to the origin (or to a given site) 
and with the mean time spent at the origin (or at a given site) [41].Unfortu- 
natelythe continuum formulation of the random walk problems is not well 
suited for treatment of these problems. Indeed, let us consider the standard 
(Gaussian-like) propagator for the random walk in d-dimensions given by 
[31] 

d 2 

G"(xi - x 2 ) = ( JL) 5 exp{-^pL} . (2.42) 

The mean time < T > spent at the origin is defined by the Laplace transform 
given by 

< T >_ Jim / iNe-»G N (0) * / . (243) 

o v ; 

Surelythis quantity is not well defined in any dimension. Use of spherical 
coordinates converts Eq.(2.43) into 

oo 

< T J dkk d ~ 3 = G(0) . (2.44) 
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More accurate lattice calculations indicate [42] that the above integral is 
actually finite for d=3 and divergent for d=l and 2. The probability 11(0) of 
returning to the origin is known to be related to G(0) as follows [41,42] : 



n(0) = 1 - m . (2.45) 

accordingly,the random walk is recurrent or transient depending upon 
n(0) being equal or lesser than one. 

Going back to our initial problem we would like to find out if the re- 
currence (or transience ) for S 1 problem can help us to establish the re- 
currence/transience for R 2 — problem. Looking at Eqs. (2. 38), (2. 39) we 
notice that, at least for zS> 1, such questions are legitimate to ask. 

Using Eq.(2.38) for A6 = we obtain, 

G (°) = v £ ^ = °- (2 - 46) 

m=— oo 

From here we conclude that the random walk on S 1 is recurrent. Let us 
consider now the full expression, Eq.(2.21) ,for the propagator in R 2 — 
plane.We obtain, 

~, oo 

0(0) = s E /*«p{-t}*.(*). 

rn=— oo q 

For integer m we get I m (z) = I- m (z). Therefore ,to evaluate G(0) we can 
use the series expansion for I m (z),m > : 



n=0 



oo / z \2n+m 



Substitution of this expression into Eq.(2.47) and use of the standard tables 
of integrals provides us with generic result: 



oo 

1=1 dxx v ~ x exp{ — — — 7x} 
J x 



2 ( P) 2 K„(2y/fo) 
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for the integral in Eq.(2.47) which ,when combined with Eqs.(2.47) and 
(2.48) leads us to the same conclusion :G(0) = oo .Using different argu- 
ments, McKean and Lyons [19] have reached the same conclusion: Brownian 
motion on once punctured plane is recur rent. This means, in particular, that 
the polymer lying in the plane cannot be entangled with another polymer 
which is perpendicular to this plane. The planarity of the problem is not too 
restrictive ,e.g. see Ref.[14].This means that the same result will remain 
correct in case of cylindrically symmetric problem. The situation changes 
drastically if,instead of one hole ,we would consider planar Brownian mo- 
tion in the presence of two holes. This problem and its generalizations is 
considered in the rest of this paper. 

III. The trice punctured sphere as hyperbolic Riemann sur- 
face 

In the previous section the task of finding the universal covering surface 
M for the once punctured plane was solved rather straightforwardly.The 
mapping z = expw (or w = ln\z\ + i(argz + 2irn)) with z G M and 
w £ M provides the desired answer. This answer was obtained, however, not 
in a systematic way but rather by a simple guessing. If two (or more) punc- 
tures are involved ,a simple guessing can hardly help. Therefore, to find the 
covering space M and the mapping function z = f{w) becomes a problem 
on its own. There are many ways to arrive at correct answer (which for the 
case of arbitrary number and arrangement of punctures may not even ex- 
ist ,e.g.see Section 4). In this work we choose the group-theoretic method 
which is directly associated with the homotopy theory used for treatment 
of diffusion on the circle S 1 discussed in Section 2. An excellent introduction 
to the homotopy (and also to the knot) theory is Ref.[43] to which we re- 
fer our readers for more comprehensive treatment(s).Some facts about the 
combinatorial group theory can be found in accessible form in Refs[44,45]. 

Let us begin with the following observation. The once punctured sphere 
S 2 — is actually simply connected(i.e.each point can be connected with any 
other so that all paths are homotopically equivalent). The twice punctured 
sphere is equivalent to R 2 — and ,therefore,the covering problem we had 
studied already in previous section. The homotopy classes in this case are 
made of paths with different winding numbers. The homotopy group is free 
infinite cyclic and abelian with just one generator a so that paths with differ- 
ent winding numbers are made of successive applications of a and a -1 . Surely, 
a m a n = a n a m for windings m and n since the group is abelian. Let us now 
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Figure 1: Homotopy of paths(with respect to some arbitrary point p) on the 
trice punctured spere 



consider the trice punctured sphere (or twice punctured plane). The homo- 
topically nonequivalent paths are depicted in Fig.l. 

The homotopy group is made of generators a and b related to the wind- 
ings around the first and the second puncture respectively. Lyons and McK- 
ean [19] and McKean and Sullivan [20] in their study of diffusion on the trice 
punctured sphere introduce yet another generator c which is associated with 
windings around the point at infinity. The homotopy group is non-abelian 
with 3 generators a, b and c subject to one relation 



The combinatorial group theory teaches us [44] ,that any group G can be 
described in terms of set of generators {aj} and relations {R{\ between some 
(or all) of group elements. The generators and relations form a presentation 
for the group G ,i-e. 



In the cases when there are no relations the groups are called free. It can 
be proven that any free group with two or more generators is non-abelian. 
It also can be proven that any group G can be considered as some subgroup 
of the free group (with sufficient number of generators). In our case we have 



abc = 1 



(3.1) 



G —< ai, ...a n | R\, Rk > • 



(3.2) 



G =< a, b, c | abc = 1 > 



(3.3) 
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This is a special case of the so called triangle groups [45] to be further 
discussed in section 7. At first sight it looks like this group is not free. This 
, however, is not the case. Indeed, let us notice first that c=(ab)~ 1 .Consider 
now the combinations ab and ba .If the relation abc = 1 would be absent 
and,instead ,we would have just free group made of two elements a and b 
,then we can construct a word W (actually ,the reduced word) via 



W = a ai b fSl ...a ar b fSr 



(3.4) 



where a\ and (5 r can be any integers while ctj and 0i are any integers 
, except zero.Evidentlythe word W encodes some specific pattern of wind- 
ings around the first and the second puncture. If,instead of elements a and b 
we would consider the combinations a = ab and b = ba , then, we can again 
construct some wordjust like in Eq.(3.4).The totality of all words based on 
generators a and b happens to be the same as that based on generators a 
and b [44]. But c is just a -1 ! Because of this,the group based on presenta- 
tion given by Eq.(3.3) is actually free group constructed of two generators, 
i.e. presentation given by Eq.(3.3) can be replaced by the equivalent presen- 
tation 



G =< a,b |> . 



(3.5) 



It can be shown[44], that such defined free non-abelian group can be realized 
with help of ,say, 2x2 matrices 



so that 



a = 



1 

x 



1 

nx 



and 



and 



1 




1 





x 
1 



nx 
1 



(3.6) 



(3.7) 



where x is some rational number. The above introduced matrices are unimod- 
ular.This means the following. Let z be some complex number. The group 
of unimodular transformations is made of transformations of the type 



az + b 

z — > : 

cz + a 

such that the determinant DetL = ab — be of the matrix 



L = 



a 
c 



b 
d 



(3.8) 



(3.9) 
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Figure 2: The fundamental triangle on H 2 



is equal to one.Clearly,transformations denned by matrices a and b are uni- 
modular. It is known [46] that the unimodular transformations are acting 
on the upper half plane( Poincare model H 2 oi the hyperbolic plane )defined 



by mapping it to itself. 

The H 2 is the universal covering space for every Riemann surface of 
genus g > 1[46, 47]. Such surface can be constructed with help of some fun- 
damental polygon D in H 2 .The appropriately chosen(see below) unimodular 
transformations will translate(tesselate) this polygon across H 2 so that the 
entire H 2 will be covered by copies of this polygon without gaps. This 
picture is an elementary extension of a much simpler picture of torus as 
a quotient of the Euclidean plane by the group of translations tesselating 
the fundamental quadrangle or of a circle S l as quotient R/Z which we had 
studied in Section 2. 

For the trice punctured sphere we can choose one puncture at infinity 
while two others at points and 1 respectively. Then,the fundamental trian- 
gle is depicted in Fig. 2. 

Naturally,the vertices of this triangle should not belong to the triangle 
itself. If the ff 2 plane is filled by such tesselated triangles, the w-plane is 
covered by an infinity of the upper and lower half-planes which are the con- 



by 



H 2 = {x + iy 



zeC\y>0} 



(3.10) 
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Figure 3: The fundamental quadrangle on H 2 .It is also used for description 
of the punctured torus, e.g.see Ref.[] and Fig. 13 below. 



formal images of the tesselated triangles. Each half-plane has three "neigh- 
bors" which are connected with it along the stretch < w < 1 and two 
rays :-oo < w < and 1< w < oo respectively.The totality of half-planes 
connected with each other in the way just indicated forms the modular sur- 
face M (one way of making Riemann surfaces [48]). The modular surface 
M is the universal covering surface for the trice punctured sphere. More 
exactly, the motion of some Brownian walker on the twice punctured plane 
around puncture(s) can be associated with the motion of its image on H 2 
.The topology of the twice punctured plane model is actually the same as 
for the once punctured torus .This fact was briefly discussed in our earlier 
work , Refs[l,14] . 

Let us discuss now some consequences of this observation. In the next 
section we are going to find the explicit form of the function which maps 
the complex plane with cuts (or better M )into H 2 . Let w = J(z) be such 
function so that w 6 M and z 6 H 2 .Consider action of such defined function 
on the fundamental triangle depicted in Fig. 2. without loss of generality ,it 
is more convenient to consider two adjacent triangles as depicted in Fig.3. 

These triangles make up the region Rc H 2 bounded by two vertical lines 
Rez=±l and two semicircles |z± ^| = ^. The left boundary Rez = — 1 is 
connected with the right boundary Rez = 1 with help of translation: b(z) = 
2 + z while the lower left boundary | z + \ | = \ is connected with the lower 
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right \z — 2 1 = 2 °y means of transformation a(z) = -Comparison with 
Eqs(3.6) and (3.8) indicates that the combination of such defined motions 
do constitute free non-abelian group composed of two generators. For the 
function J(z) we have to require 

J(a(z)) = J(z) (3.11a) 

and also 



J(b(z)) = J{z) (3.11b) 

in order to achieve the tesselation of H 2 without gaps. It could be shown as 
well that 



J(ab{z)) = J(ba(z)) = J(z) (3.12) 

Such equations are the defining relations for the automorphic function J(z). 
At this point, several questions arise. First, are the above transformations 
a{z) and b(z) unique ?Second, based on these transformations alone can one 
find the explicit form of the function J(z) by using the periodicity properties 
shown above ? The answer to the first question is negative. Already Eqs.(3.6) 
and (3.7) indicate that this is not the case. At the same time, for given set of 
a and b generators it is possible to find the automorphic function J(z) which 
will possess the required transformation properties. Non uniqueness of the 
transformations leads to the non uniqueness of the fundamental polygon 
D on which these generators act. Recent studies had shown [49] that some 
choice of generators of the free group G(/j,),e.g. 

b \ Z \ = 2+ A (3.13) 

where fii = 3f;or ^2 = 0.0533 + 1.9i ,etc also describe the homotopy group 
of the punctured torus. For the case fj,i = 3i a simple connected domain 
De H 2 is depicted in Fig.4. 

Action of generators a(z) and b{z) on the fundamental polygon does 
not lead to the tesselation of the entire H 2 .Instead ,the above free group 
tesselates some subdomain A of H 2 so that the quotient A/G is still the 
punctured torus. For fi = ^2 the domain A is depicted in Fig.5 and has, 
by now familiar, fractal-like shape.We shall not go into details of such much 
more elaborate treatment of the punctured torus .For less exotic situations 
, e.g. that depicted in Fig. 3, the explicit form of the function J(z) can be found 
and this task is accomplished in the next section. 
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Figure 4: A simply connected domain D of H 2 whose quotient D/G(//i) 
produces the punctured torus [] 




Figure 5: A simple connected domain D of H 2 whose quotient D/G(//2) 
produces puncured torus as well.f]- 
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IV. Uniformization of the n-punctured sphere 



On the sphere S 2 = C U {00} let us consider a region $7 obtained by 
removing a finite number of points, i.e.Q = S 2 \{pi, ...,p n } with n> 3 .Ac- 
cording to the uniformization theorem there must be some meromorphic 
function A which provides the universal covering map of £1 by H 2 defined by 
Eq.(3.10). Suppose J is some function which maps H 2 into Q ,then if g is 
some function which maps H 2 to H 2 (e.g. see Eq.(3.8)), then J = \og , where 
o denotes the functional composition. Surely, such defined function J coin- 
cides with that which was introduced earlier, e.g.see Eq.s (3. 11), (3. 12). The 
function J is locally invertible so that a closed path in f2 has as its image a 
closed path in H 2 if and only if it is homotopic to zero. A closed path in Q 
with ,say ,winding number +1 with respect to p k and with respect to pj, 
j / k , lifts to a path in H 2 with the point T k z ,where is related to one 
of the homotopy generators ,e.g. a (or b ) , discussed in the previous section. 
Evidently, by analogy with Eq.(3.1), we can write in this more general case 



Let us now connect points piand P2 by a simple (without crossings) arc q± 
and let us proceed in the same fashion with P2 and p% , ...,p n and p\. Then,fi 
is being decomposed into two simply connected regions G\ and Upon 
mapping into H 2 these regions will go into subregions of H 2 whose boundary 
meets the real line {y = 0) at 2n-2 points. Fig.3 represents just an example 
of this sort. In this case we have 3 points -1,0 and +1 of the fundamental 
region R located on the real line while the fourth point is located at ±00 of 
real line. In order to find such mapping the following theorem is helpful[50]. 

Theorem.4.1. Suppose z = J{w) where w £ H 2 is universal covering 
map of {l,then,{w,z}is meromorphic function of z given by 



with w = w(z) = J 1 (z) and w' = j]j,etc ,and constants are subject to 
relations: 



T x T 2 ...T n = 1 



(4.1) 




(4.2) 



where the Schwarzian derivative {w, z) is defined by 




(4.3) 



n 



m k = 0, 



k=i 
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n 

J2( 2 m kPk + 1) = 0, (4.4) 
fc=i 



k=l 

Eq.4.2 is the third order nonlinear differential equation whose solution may 
or may not be easy to find. However, already Poincare' had discovered [29] 
that ,instead of struggling with solution of Eq.(4.2) ,one can consider much 
simpler second order linear differential equation of the type 

y" + ±{ w ,z}y = 0. (4.5) 

In order to solve Eq.(4.5) the r.h.s.of Eq.(4.2) must be used. Then, the ob- 
tained equation is of Fuchsian type [30] and can be solved, provided that the 
constants m k (the accessory parameters) are known. Eq.(4.5) produces two 
linearly independent solutions y\ and y2-The most spectacular outcome of 
all this can be summarized as follows: 

Theorem 4.2 . If two independent solutions of Eq.(4-5) are known, then 
the mapping function w = w(z) can be determined from solutions of one of 
the following equations: 

yi = w (w/) ^ or i/2= («/) 5 . (4-6) 

To find the accessory parameters m k is not an easy task in general. 
Recent results of conformal and string theory provide some new results in 
this direction. Before discussing these results, let us illustrate how the results 
just described work in the case of twice and trice punctured sphere. In the 
first case we may associate p\ with z = so that the combined use of 
Eqs(4.2)-(4,5) produces the following differential equation: 

y" + = o • (4-7) 

Looking for solution of Eq.(4.7) in the form y = z a we can easily obtain 
Vi,2 = \fz .Using Eq.(4.6) we obtain as well 

v^=(7^ 7 M) _1 (4.8) 
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This leads us immediately to 

w = In \z\ + const (4-9) 

The obtained result is incomplete agreement with earlier obtained,Eq.(2.16). 

In the second case we may associate p\ with z = 0, P2 with z = 1 and 
Pz with z = oo. The combined use of Eqs(4.2)-(4.5) then produces 

v" + {\+ At 1 TT2 +7-- 71 1 TT >y = - 
z 2 4(z-l) 2 4z 4(z-l) 

This equation can be reduced to that for the hypergeometric function if we 
use the following substitution: 

y(z) = y/z(z-l)f(z). (4.11) 

Then, for the function f(z) we obtain the following hypergeometric equation: 

z{z-l)f" + {l-2z)f -\f = Q. (4.12) 

Let us recall that the hypergeometric function F obeys in general equation 
of the type 

z(z - 1)F" + [c - (a + b + l)z]F' - abF = (4.13) 
which has solution in terms of the following series expansion 

„, , . a- b a • (a + 1) • b • (b + 1) 9 . , ,. 

F(a, t , q ,) = l + — z+ l {c ; i} ,\, 2 > * + -.. (4.14) 

For a special values of a, 6 and c and also z this expansion coincides exactly 
with the series expansion for the complete elliptic integral K defined by 

7T 

2 



K = / ^ / „ ■ 2 • (415) 

J vl-K sin x 
v 

For this to happen,we have to require a = b = \,c = 1 and z = k 2 in 
Eq. (4. 13). This then produces [48] 

K=|F(i i l;fe 2 ) (4.16) 
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Using these results we conclude that K is solution of Eq.(4.12) (since the 
values of constants a, b and c produce Eq.(4.12)) in which we have to replace 
z by k 2 . Another independent solution of Eq.(4.12) can be obtained if we 
replace k 2 in Eq.(4.16) by (k 1 ) 2 = 1 — A; 2 . From the theory of elliptic functions 
it is known that (A/) 2 should not touch the real axis which is cut from z=0 
to z=-oo .This circumstance implies then that k 2 should not touch the real 
axis cut from z=l to z=+cx) .Using Eq.(4.6) we obtain now 

where K (z) is the same as K in Eq.(4.16) with k — > k' and k! being replaced 
by z .Strictly speaking, solutions of Eq.(4.13) are determined with accuracy 
up to some constant c.Therefore,the above result should be amended by 

w{z) = c^ = J~ 1 (z) (4.18) 

The choice of constant c is determined by the requirement that the above 
mapping of z-plane (with cuts from 1 to oo and from -oo to 0)is mapping 
onto the fundamental region of the Poincare H 2 plane as depicted in Fig.s 
2 and 3. If we choose c=i ,then the upper z half-plane is mapped into the 
right half of the fundamental region D depicted in Fig. 3 while the lower z 
half plane is mapped into the left half of D [51]. The inverse of w(z) is given 
by the function J(w) defined below: 

oo / i i 2n \ 8 

z = J(w) = 16, nirT^rj ( 4 - 19 ) 

n=l ^ ' 

with q = exp{iirw}. These results were briefly mentioned in our earlier 
work,Ref [14],e.g.see Appendix of Ref.[14]. Generalization of the obtained 
results to the case of n punctures (n> 3)is nontrivial (this is so because the 
system of equations, Eq. (4.4) can be easily solved only for n=l and 2, pro- 
vided that the third point is located at oo). Therefore, many methods can be 
used, in principle, with variable degree of success. In this section we would 
like to discuss only methods which originate in string and conformal field 
theories. In section 7 we shall discuss in some detail other possibilities. 

In general, the accessory parameters mk = mk(pi, ■■■,Pk) could be some 
very complicated and, in general, unknown functions of puncture locations. Some 
progress had been recently made based on the following key observation 
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made by Zograf and Takhtadjian[52]. Consider the Schwarzian derivative, 
Eq.(4.3), rewritten in the following way 



{w,z} = -- 



d z ln 



dw 



d2 , 



dw 



(4.20) 



This result coincides with that given by Eq.(4.3) as can be seen by direct 
calculation.At the same time,let <p(z) = In |^|,then we can write instead of 
Eq.(4.20) the following result 



d 2 1 (dp> 



(4.21) 



or,in view of Eq.(4.2), 



dz 2 



r 



dip 
dz 



1 n 



+ 



fc=i 



Pk 



(4.22) 



The function ip satisfies the Liouville equation 



a 2 

dzdz 



-exp<p(z,z) 



(4.23) 



which describes the surfaces of constant negative curvature k = — l.The 
Poincare H 2 model ,Eq.(3.10), equipped with such type of metric becomes a 
model of the hyperbolic (or Lobachevski ) plane. The Liouville equation can 
be obtained as equation of motion coming from the variation of the Liouville 
action functional S[(p] given by 



S[<p] 



d 2 z 



+ counterterms 



(4.25) 



where A is determined by A = fi\ Li™!} 1 {\z — Pi\ < r} U {\z\ > r — ► + 
and the explicit analytic form of the counterterms is not important for what 
will follow next. By swithching from the Lagrangian to the Hamiltonian 
formalism we obtain , using Eq.(4.25) , the energy-momentum (the stess- 
energy) tensor T zz given by 



d 2 

dz 



2<P 



dp 

dz 



= T 



(4.26) 
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This result can now be combined with Eq.(4.22). This leaels,in view of 
eq.s(4.5), (4.21) and (4.22) to the result: 

y" + l -Ty = 0. (4.27) 

This equation represents the classical limit of the corresponding quantum/statistical 
mechanical equation considered in the famous paper by Belavin , Polyakov 
and Zamolodchikov (section 5 of Ref.[53]) in connection with the develop- 
ment of the conformal field theory. Systematic corrections to these classical 
results ,Eq.s(4.20)-(4.27), can be achieved through Feynman's way of doing 
quantum mechanics. That is one formaly defines the path integral 

<n>= j Z%]exp{-^M} (4.28) 

where Q = S 2 \{pi, ...,p n } as before, C(f2) is the appropriatly chosen domain 
of functional integration and h is some control parameter which is fixed at 
the end of calculations [54]. It is not our purpose here to discuss the impli- 
cations of such path integral formulation of the uniformization problem. The 
details could be found in the references already provided. Here we only are 
concerned with calculating the accessory parameters m^. These parameters 
can be obtained ,in principle, by means of a simple looking formula[54,55]: 

m k = (4.29) 

The simplicity of this formula is somewhat mileading since it was obtained 
with help of a saddle point approximation .And even within this approxima- 
tion actual computation can be made only for some special cases,e.g. when 
two punctures are very close to each other while the rest are far away. In 
particular, for n<4 ,one obtains [55] 

m k = -7u — " V' Pk ~^Pn, (4.30) 

2(p fc -Pn) 

while for n> 4 and one of the punctures is being located at zero ,one obtains 



1_ vr_ 2 

2p 2p(ln\p\Y 



m = -— + - - ; 2 , p — >0. (4.31) 



It is also possible to obtain the exact results for some very special symmet- 
ric arrangement of the punctures. Whence, although, in principle, solution 
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of Eq.(4.5) solves the uniformizaton (or mapping) problem completely, in 
practice, even if all would be known, solution of Eq.(4.5) could be so 
complicated that its practical use in path integral calculations similar to 
that described in section 2 becomes nonpractical.Nevertheless,our efforts 
are not completely in vain. For example, we have learned already from 
example of the trice punctured sphere that solution of Eq.(4.5) not only 
provides the desired map but also gives the exact shape of the fundamen- 
tal region in H 2 plane , e.g. see Fig.s 2 and 3. In order to use other meth- 
ods,e.g.see section 7,we also need the precise shape of the fundamental re- 
gion or,equivalently,the explicit form of the unimodular transformation, Eq. 
(3.8). For the fundamental domain depicted in Fig. 3 the corresponding uni- 
modular transformation is known to be [51] 

1 2 
1 

in accord with earlier obtained Eq.(3.6) .Analogously, if the points are ar- 
ranged as n-th roots of unity , it can be shown [50], that 

m k = - l - V - k (4.32) 

and the corresponding unimodular transformation can be found explicitly 
as well. Fortunately, this is not needed as results of section 7 indicate. The 
most important for us is the fact that the results obtained in this section for 
the trice punctured sphere could be used for unformization of the Riemann 
surface of any genus .In this work we would like only to provide an outline 
of these far reaching results, in section 7. In the meantime, in order to keep 
things in proper perspective,we would like now to discuss some topics which 
are directly related to just covered. 



a = 




1 



and b = 



V. Connections with the Riemann-Hilbert Problem and 

Knizhnik-Zamolodchikov equations 

Eq.(4.5) with {w,z} given by Eq.(4.2) is a typical case of Fuchsian dif- 
ferential equation []. The most general form of equations of Fuchsian type is 
known to be given by 

y<P> + qi{x)y^ + ■■■ + q p (x)y = (5.1) 

where the coefficients qi{x) near some singular point a are given by 

qi(x) = -^L, i = l,...,p (5.2) 
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where functions r±(x), ...r p (x) are holomorphic at a. The above p-th order 
linear differential equation has, in general, p independent solutions .To find 
these solutions it is convenient to replace the above higher order differential 
equation with the equivalent system of the first order differential equations. 
This can be accomplished through the following changes of variables : 

V = z 1 , 



< \ d y 

{x ~ a) Tx 



(x — a) 



p-i 



y 



dx p 1 x 

Using such defined z variables and the fact that Eq.(5.1) is equivalent to the 
system of equations given below 

dy 1 2 
-dx- = V 



dx 



(5.3) 



djf 
dx 



= -q P (x)y 1 qi{x)y p 



it is a simple matter to bring Eq.(5.1) to the canonical form 



Tx = A{x)z ' 



(5.4) 



where the vector z is given by 



(5.5) 
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and the matrix A(x) is given by 



A(x) 









1 
1 







1 

2 





1 



x — a 



r v— 1 







p — 1 — r\ 



(5.6) 



It is clear that the above arguments can be extended to the case when 
, instead of one singularity located at a we would have n singularities located 
at di , i=l-n. Then, instead of Eq.(5.4),we would have equation which looks 
like Eq.(5.4) but with matrix A(x) being replaced by 



A(x) = £ 



i=l 



x — a,- 



-Bi(x) 



(5.7) 



with matrix Bi{x) being nonsingular at above points. For example, the hy- 
pergeometric equation, Eq.(4.13), can be presented in Fuchsian form given 
below 



dz 
dx 





-ab 







1 



+ 



x 





c 



1 

(a + b) 



x-1 



(5.8) 



while the Bessel equation , Eq.(2.34), which upon trivial rescaling can be 
brought to the standard form 



d 2 y 1 dy 

dx 2 x dx 







can be presented in the Fuchsian form of the following type 

1 



dz 
dx 



v 2 - x 2 



where the column vector z is given by 



z = 



y 

dy 

J dx 



(5.9) 



(5.10) 



(5.11) 



Eq.(4.7) is also of Fuchsian type . It can be brought to from of Eq.(5.10) 
with matrix A{x) given by 



A(x) = - 
x 



1 
-1 1 



(5.12) 
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with vector z being the same as in Eq.(5.11). 

Consider now Eq.(5.4) from different perspective. First , assume that 
there is some matrix function T(x) such that we can relate vector z to 
another vector Z via 



Z = T(x)z. (5.13) 

This then allows us to write 



or,equivalently, 



dZ _ dT ^dz 

dx dx dx 



^ = r-^-r-^z. (5.14) 

dx dx dx 



Substitution of this result back into Eq.(5.4) and elimination of Z with help 
of Eq.(5.13) produces 

dz 

-7- = A'(x)z (5.15) 
dx 

with matrix A'(x) being given by 

7T-1 

A'(x) = TAT' 1 + ^-T' 1 . (5.16) 

dx 

At this point the reader can easily recognize the gauge transformation char- 
acteristic of the nonabelian gauge field theories [56]. Unlike the field-theoretic 
case, our field is nonrandom and,therefore, we cannot use field-theoretic 
methods which always involve the averages over randomly fluctuating fields. Such 
averages make sense , nevertheless, since in physical reality the punctures(poles) 
are not fixed but can move(see, however, section 7 for further details).. One 
may think of Eq.(5.4) (or (5.15)) as some sort of equation for the massless 
Dirac particle(e.g.neutrino) placed in time-dependent (or space-dependent 
but not both!)nonabelian gauge field [56]. The question arises: can one con- 
sider instead the massive Dirac particle in such gauge field ? The answer 
turns out to be negative as we are going to demonstrate momentarily this 
purpose instead of Eq.(5.7) we would like to consider more general matrix 
given by 



n 

A{x) =J^^— Bi{x)+B{x) (5.17) 



-i x 
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where the matrix B{x) is not singular for all x . Let all a\ / oo and require 
our system of equations ,Eq.(5.7) to be nonsingular at oo as well. This puts 
some constraint on both Bi(x) and S(x).By replacing x with r = x _1 in Eq.s 
(5. 7), (5. 17) a simple analysis shows that B(x) = and, in addition, 

n 

J2Bi(x) = (5.18) 

i=i 

Accordingly, we cannot use massive Dirac particle interpretation of the ob- 
tained results. At the same time, let x £ O where £1 was defined at the 
begining of section 4. Let ,as before,-?/ 2 be the universal covering space of Q 
(that is we are considering at least three punctures). By analogy with sec- 
tion 3 ,let G be the group of motions in the covering space. Let now a and 
b £ G ,then if z{x) is solution of the system Eq.(5.7) and z{x) is solution in 
the covering space,then as before, z(ax),is also a solution to Eq.(5.7) in the 
covering space .Moreover, in view of the linearity of Eq.(5.4) ,such solution 
is defined only with accuracy up to some constant matrix x( a )- Let us in- 
troduce notations:z(<7x) = z o a and z{ax)x{ a ) = {zoo) x( a )- Using these 
results,we obtain as well 

z o 5 = [(z o a) x( a )] o 5 = (z o a o 5) x(cr) 

This naturally implies 

(z o 5) x(S) = (zoao5) x(cr)x((5) = (z o (aS)) x{v5) 

that is 

xM=xWxW (5.19) 

which is a monodromy group condition analogous to that we have already 
encountered in connection with Eq.(2.7). In case of the trice punctured 
sphere, using Eq.(3.1), we obtain 

X(abc) = X (a)x(b)x(c) = x(l) = 1 (5.20) 
Evidentlyin case of more than 3 punctures we should have 

n 

IJx(a<) = l (5.21) 

i=i 

where is the generator of motion around i-th loop (e.g. see Fig. 1). The 
classical Riemann-Hilbert (R-H) problem can be formulated now as follows. 
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R-H problem: given location of the points p\, ...,p n G S 2 and matrices 
xifli), find a linear differential equation of the type given by Eqs(5.4)--(5.7) 
whose monodromy group coincides with the group associated with matrices 

X(a»)- 

If we allow points p\, ...,p n to move, then one may pose 
The related problem: find conditions under which for the fixed mon- 
odromy group (and this is very natural since topologically nothing changes) 
solutions of Eq.(5.4) and the matrix A(x) will depend continuously on the 
set pi, ...,p n . 

More precisely, let z = z(xq;x) be solution of the system of Eqs(5.4) 
,where xo = {pio, ■■■,Pno} -Then, the necessary and sufficient condition for 
the matrix A(x) as function of parameters pi,...,p n and xq to be associ- 
ated with the fixed monodromy group depends upon the solvability of the 
following set of (consistency) equations 

^i = [A j ,A i }(— 

OPi Pi - Pj PiO - Pi 

^i = -Y^[A i ,Aj]^—, (5.22) 

dpi 4rr Pi- pj 

dAi v-^ r . 1 . 

— = y\Ai,Aj\ ,i = l, ...,n . 

dx f-' x -pj 

The system of equations given above is known in the literature as Schlesinger's 
equations [57]. Study of solutions of these equations is directly associated 
with study of the exactly integrable classical and quantum systems [58] , 
conformal field theory[59] and Einsteinian 3+1 gravity[60] and , therefore, is 
not going to be discussed further in this section. 

Nevertheless,to make our presentation self-contained, we would like to 
discuss briefly connections between the system of Eqs.(5.7) and that known 
in the literature as Knizhnik-Zamolodchikov (K-Z) equations [61]. The sys- 
tem of equations (5.4) can be formally solved using method of Lappo- 
Danilevsky [62]. Let the matrix of solutions of Eq.(5.4) be just a unit matrix 
I for x = Xo,then , for x ^ xq we can formally write solution to Eq.(5.4) in 
the following form 

oo l,...,n 

z(x) =I + Y^ 2 L xo(Pjl,~;Pjv I X ) A h ■ ■ ■ A J V ( 5 - 23 ) 

"=1 
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where the functions L XQ (...) are defined recursively as follows 



L Xl A,> n ■<•) / (5-24) 
y Pji 



Xo 



and , accordingly, 



LrJ.,P n ,> h I / PA, - l|y) rfy- (5-25) 



Using these results, let us consider now the following system of equations 



9 £ = h± J^rfi-n <•„)• (5.26) 



Here /i is some prescribed constant (control parameter) and the symmetric 
matrix is also assumed to be given (actuallythe matrix is not the 
matrix but rather the tensor product of two matrices [61] but for the mo- 
ment we shall ignore this fact ).The above system of equations is known 
in the literature as K-Z system of equations [61]. It has numerous physical 
applications , e.g. see Ref.[14] , which are not going to be discussed here. What 
will be important for us is its connection with Eqs(5.4) and (5. 22). To es- 
tablish this connection several steps are needed. First,using the property 
of symmetry of the matrix t{j it is sraightforward to obtain the following 
auxiliary equations 

E£ = 0and I>£= E /(*i.->*n) (5.27) 

i=l * i=l * l<i<j<n 

Eq.s(5.27) reduce the number of independent variables in K-Z equation to 
n — 2 .In particular, solution f(x±, X2, x 3 ) of K-Z Eq.(5.26) actually depends 
only upon one variable which we shall denote as x.More specifically, if we 
look for solution of Eq.(5.26) in the form 

f(xi, x 2 , x 3 ) = (x 3 - Xl ) h ^ 2+t23+t ^g(x) (5.28) 

where 

x 2 - X\ 



x = 



rr 3 - Xl 
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then,substitution of such an ansatz into Eq.(5.26) produces the following 
equation for the function g(x): 

d f = h(^ + ^-) 9 { X ). (5.29) 
ax x x — 1 

In the original paper by Knizhnik and Zamolodchikov[61] it is explained 
why ti2and £23 are actually 2x2 matrices (not just usual numbers which are 
components of the matrix). Comparison with Eq.(5.8) then suggests imme- 
diately that we are dealing with the hypergeometric equation and ,indeed, 
Knizhnik and Zamolodchikov had obtained two solutions for g{x) in terms of 
hypergeometric functions .Their solution is not the most general , however. 
A complete solution was found by Drinfeld[63] and ,more explicitly, by Le 
and Murakami [64]. It is noteworthy to describe it since it is related to quan- 
tum groups and Vassiliev invariants for knots and links [65-67]. To describe 
this solution, let us consider instead of Eq.(5.29) the following generalization 

^ = — ^=(- + AM-) MO) 
dx 2ttV^T v x x- J y ' 

where the prefactor 2 ^ = is chosen for convenience only (since ^ is just 
the usual Planck constant in quantum mechanics)and the factors A and B 
are some non-commuting variables, e.g. matrices. Let g\{x) and g2{x) be two 
independent solutions of Eq. (5. 30). Then, as it was shown by Drinfeld [63], 
gi(x) = g2(x)&(A, B) where the function <E> of two noncommuting variables 
is known in the literature as Drinfeld associator .The explicit form of <I> is 
hard to find in closed form .If h is small parameter, then , after some tedious 
calculations [64] which involve use of Lappo-Danilevskii-type expansions ,one 
finally obtains 



(2^1) 21 (27T y^iy 

(5.31) 



*(A,B) = 1 - - ^> [A,B]h? + —±->([[A,B],B] - [A, [A,B]])h 3 + ... 



Here q(x) is Riemann zeta function. To obtain solutions gi(x) and 52 (x) is 
rather easy in the neigbourhoods of and l.To this purpose it is sufficient 
to look for solution of g\ in the form 

gi{x) = P(x)x hA (5.32) 

and, analogously, 

g 2 {x) = Q{1 - x)(l - x) hB (5.33) 
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Assuming that functions P(x) and Q(l-x) can be represented in terms of 
power series expansions it is rather straigthforward to find the successive 
coefficients. For values of x not necessarily close to x=0 and x=l again one 
should be looking for solution using Lappo-Danilevsky method. Altschuler 
and Freidel had demostrated [66] how Drinfeld associator $ can be used 
for calculation of the regularized Vassiliev-Kontsevich invariant for knots 
and links. It is not our purpose to reproduce these calculations in this pa- 
per. Instead,we would like to mention yet another physical application of the 
obtained results. 

To this purpose, let us consider again system of equations Eq.(5.4) with 
matrix A(x) being given by Eq.(5.7).Let us first choose matrices Bj to be 
constants. Then ,the gauge transformation ,Eq.(5.16), is just a similarity 
transformation so that we can always select such F's that the resulting 
matrix be diagonal. Without loss of generality ,we shall assume that this 
is the case from now on. Let us then look for a solution of Eq.(5.4) in 
the vicinity of one of the points, say, aj. We can always choose the system 
of coordinates located at this point. Then, in the vicinity of this point we 
essentially have to solve the following equation 

— = -*-z (5.34) 

ax x 

The solution of this equation is obtained at once and is given by z = x Aj 
=exp{A,- In a;}. Now, if x is in the complex plane (or on S 2 ) use of Eq.(2.16) 
produces , after one rotation around x=0, the following solution z o a = 
cxp{2niAj}z = x( a ) z - If now the matrix A(x) is x-dependent,then one can 
look for solution of Eq.(5.4) in the form 

n 

z {x) = ®i( x )( x ~ a i) Aj ( 5 - 35 ) 

3=1 

where Aj is the value of A(x) at x=dj (surely ,upon proper diagonalization) 
and the functions <&j(x) are holomorphic at x=cij .Instead of S 2 ,let us con- 
sider a disc D 2 first. Let xo be the location of singular point(s) in D 2 and let 
71, 7 n be a set of loops which go around the singular points while { x( c7 i)} 
be the respective monodromy matrix set. Now we cover D 2 by domains Uj 
so that each domain Uj has only one singularity aj .Using Eq.(5.35) in each 
domain we have zj = &j(t)t Aj (where t is the local coordinate in U^-th do- 
main) so that at the intersection of Ui and Uj domains we should have the 
consistency condition 

$T 1 $ j+1 = t A it- A '+ 1 | Uj n U j+ i (5.36) 
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or ,more generally, if one considers S 2 , 

= G(x)^.(x) (5.37) 

with G(x) being some known function . This type of boundary value prob- 
lem is very well known in the theory of elasticity [68] and was thoroughly 
studied since it was first formulated by Hilbert in 1900 (21st Hilbert's prob- 
lem) [69]. Whence, the use of the non -Abelian gauge field theories has a 
long history . Birkhoff [70] had shown that solution of the Hilbert problem 
,Eq(5.37) is equivalent to the solution of the Fuchsian system of equations , 
Eq.(5.4). More readable account of this proof is contained in Ref.[71] .Use 
of Riemann-Hilbert problem in the the theory of exactly integrable equa- 
tions is well documented, e.g. in Ref [72]. Accordingly, it is not our intention 
to discuss these topics in this paper. Instead we would like to discuss topics 
related to random walks on groups in the following section. 

VI. Random walks on free groups: an introduction 

By now, it should become clear that even for 3 punctures in the plane 
the path integral approach developed in section 2 becomes nonpractical.In 
current literature the situation with 3 or more punctures is treated with 
help of an auxiliary abelian or nonabelian gauge fields as described in our 
Phys. Reports , Ref [14]. Such approach though plausible but mathematically 
less rigorous. To realize the difficulty , we would like to mention the fol- 
lowing. In Ref. [7] the trice punctured sphere (or twice punctured plane) 
problem was solved with help of path integral methods. However,the major 
question about the recurrence/transience was not even addressed. More- 
over,the method used in Ref. [7] does not admit generalization to the multi- 
ple punctured case . Accordingly, it is necessary to develop other methods of 
dealing with these problems. In this section we shall discuss an alternative 
approach which originates from the fundamental work of Kesten[73] . 

In section 3 we had discussed presentation for the group G,e.g.see Eq.(3.2) 
(with all Ri = 0) and associated with it word W ,Eq. (3. 4). Following Kesten, 
let us consider the random walk on G in which every step consist of right 
multiplication by a^ (or its inverse) each with probability pj. Surely, one can 
complicate matters by assigning different probabilities for direct and inverse 
elements but we shall not complicate our presentation at this time. For the 
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free group of n elements we should require at each step 

n 

2^> = 1 (6.1) 

i=i 

The random walk thus defined determines a Markov chain whose possible 
states are elements of G. The transition probability(to be precisely defined 
below) from some W\ G G to another W 2 G G is given by the probability that 
W 2 is reached in one step from Wi . So far all this looks very detached from 
the previous discussions. To make a connection, we shall , following Ref.[3] 
, introduce some additional concepts. Let X be a metric space in which the 
distance (the geodesic distance) d(x,y) between points x and y which belong 
to X is formally determined by 

d(x,y) = \x-y\ (6.2) 

In particular, if W\, W 2 G G ,then 

\Wi - W 2 \ = mm{m \W^W 2 = a ia2 ...a m }. (6.3) 

That is the distance between words is determined my the minimal word. For 
example, let W\ = a^a^ 1 ^ .Surely ,this word is not minimal since the 
combination a 2 a^ 1= 1 can be easily removed. As simple as it is, the general 
problem of comparing different words , known as word problem [44] , cannot 
be solved in general [74]. This does not mean that the problem cannot be 
solved for a specific group,for example, for the free group. More on this sub- 
ject is discussed in the classical monograph, Ref. [44]. More important for us is 
the fact that such defined distance,Eq.(4.3), possesses the property of trans- 
lational invariancejust like in the case of Eq. (2. 3). Indeed, let us consider 
\jWi — 7W2I .According to Eq.(6.3) we have 

I7W1 - 7W2I = (7W4)" 1 7W2 = WfS'SWb = W^W 2 . (6.4) 

It is possible to introduce not only the distance but even the curvature of the 
group, etc. For instance, following Milnor[75], for a finitely generated group, 
i.e. group of finite number of generators, e.g. n , one can define the growth 
function 7(5) which for each positive integer s determines the number of 
words with length < s. Then, it can be shown , that for the free group of n 
generators the growth function is given by 

7 (s) = (n(2n-l) s -l)/(n-l) (6.5) 
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and this result is directly related to the curvature of the underlying manifold 
whose fundamental group consist of n elements. 

It should be noted, that the distance defined above is not the only possi- 
bility in groups. For instance, if we fix some element, say unity, as a reference 
point, then ,one can introduce another distance between x and y ,also known 
in the literature as Gromov product [3,76,77], e.g. 

dc(x,y) = ^{\x\ + \y\-\xy' 1 \) (6.6) 

so that, by construction, dc(x, x)=\x\ > which makes perfect sence. Using 
this result one can go further by providing the following 

Definition 6.1. Let 5 be non-negative real number. The metric space X is 
said to be 5— hyperbolic if 

da{x,y) > mm(d G {x,z),d G (y,z)) - 6 (6.7) 

for every x,y and z£ X. Moreover, the metric space is hyperbolic (in the 
sence discussed in sections 3 and 4) if there is a real number 5 so that X is 
S— hyperbolic. 

It can be demonstrated[3,76,77] that: 

a) every graph is 0— hyperbolic; 

b) every free group of finite rank is hyperbolic. 

Moreover, following Ref.[77], the 5— ultrametric space X can now be 
defined via 

\x — y\ < max(|x — z\ , \y — z\) + 5 (6-8) 

for all x,y,z £ X .Using such definition, the following theorem can be 
proven [77] 

Theorem 6.2. Every metric space X which satisfy the above inequality 
is actually 5— hyperbolic. 

This means that the ultrametric space is a special case of hyper- 
bolic space. This fact could be used for solving some problems related to 
spin glasses and other systems which require use of replicas [4]. 

Going back to our problem of symmetric random walk on free groups we 
need to introduce two probabilities : 

m (n) _ Probability of returning to the unit element e at n-th step given 
that one initially starts at e. 

r (n) =p ro bability of returning for the first time to e at n-th step at 
n-th step given that one initially starts at e. 
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With such defined probabilities the generating functions m(x)and r(x) 
can be now formally introduced as follows: 



m{x) = ^ m^x 11 , provided that = 1 (6.9) 

n=0 

and 

oo 

r(x) =Y / r {n) x n . (6.10) 



n=l 



As it was demonstrated by Feller [78] ,e.g.see page 311 of Feller's book. 

„,(*) = T -1 R . (6.11) 

Thus,if r(x) is known, m(x) is known as well. To determine r(x) one needs 
to determine explicitly. Elegant solution of this combinatorial problem 
is presented in Kesten's paper [73], e.g see pages 348,349 of Ref. [73] and take 
into account similar derivations in Feller's book, Ref [78], with the following 
result: 



i 



r(x) = *-(" 2 -(2"-l)* 2 )*. (6 . 12) 
w 2n-l 

Use of Eq.(6.11) provides us with the following value of m(x) 

2n — 1 

m(x) = j- (6.13) 

n — 1 + (n 2 — (2n — l)x 2 )2 

with n being the number of generators in the free group G. 

Now, we would like to use all this information in order to discuss the 
random walk on Bethe lattice ,Fig.6, which is just the universal covering 
space of figure eight (in the case of free nonabelian group of two generators) 
or, more generally, the covering space of the union of n circles [74] (in case of 
n generators). Such problem was solved by conventional methods in Ref. [79] 
and we would like to borrow a couple results from this reference for the sake 
of comparison. In particular, let Pt{l |m)be the probability that the walk 
originated at position m will end up at position I after t steps. In case of the 
Bethe lattice the position m (or I) is determined with respect to the origin, 
e.g. see Fig. 6 ,so that m means the number of links making a path between 
the origin and point to, etc .The Markov chain equation is given by 

P t+1 (l\m)=Y,l{lX)Pt(l'\m) (6.14) 
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Figure 6: Figure 8 and its universal covering space 



with 7(Z, I') being determined by 

7 (J, /') = ( , (1 " + i - 1 (6.15) 

L , i = u 

The probability 11(0) of returning to the origin , defined in Eq.(2.45), is 
calculated to be 

n(0) = ^ (6.16) 

and is surely less than one for n>2. Hence, for n> 2 the walk is transient and 
this result is in accord with results obtained by Lyons and McKean[19] and 
McKean and Sullivan[20] who used different methods . For n=l the walk 
is recurrent, since n(0) = 1, and this result is in accord with earlier ob- 
tained in section 2 for 2 punctures one of which being located at infinity. We 
would like now to reproduce this result using the theory of random walks 
on groups. This is needed for various reasons discussed in the Introduction 
, below, and to be discussed in the next section. To begin, following Kaste- 
leyn[23] ,we consider an analogue of Eq.(6.14) for group elements (letters) 
ai,...,a n of some free group G . Let Pt(W) be the probability of creat- 
ing the word W of t letters starting from e. The Markov chain equation 
analogous to Eq.(6.14) can now be written as 

Pt{W) = Pt-i(W- x W')p(W) ,t > 1. (6.17) 
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with p(W) being denned as stepping probability which is the probability of 
choosing a particular branch at each vertex of the tree. Looking at Fig. 6 ,we 
can easily notice that it is exactly the same probability as for crossing the 
vertex of figure 8 (in case of two generator group). Surely, one can complicate 
matters by assigning different weights to different choices of crossing the 
vertex as it is done ,for example, in the case of the Ising model. Such 
possibility immediately connects the walks on groups with various statistical 
mechanical models as discussed ,for example, in Ref. [80]. We shall in this 
work restrict our discussion to the simplest possible case ,e.g. p(W) = 
l/2n (motivated by Fig. 6). As in the case of usual random walks [42], it is 
convenient to introduce the generating function P(W, z) as follows 

oo 

P(W, z) = p t{W)z t (6.18) 
t=o 

This generating function obeys an equation of the following type 

P(W',z)- 8 w , je = z Pt-i(W- 1 W',z)p(W) (6.19) 
weG 

so that the probability II(W) that a walk starting at e will ever visit element 
W (excluding the start as visit to e) is given by 

U(W) = km (6.20) 

v ; *-i P(e,z) 

The crucial moment now is to notice that , actually, P(e, z) = m(z) where 
m(z) is given by Eq.(6.9).At this moment Kasteleyn claims that this obser- 
vation allows one to determine P(W,z) with help of Eq.(6.19). Since the 
details of the derivation are not provided in Kasteleyn's paper, we restore 
these calculations in this work. To this purpose, let us assume that solution 
could be sought in the form 

P(W,z) = X m P(e,z) (6.21) 

with \W\ being defined by Eq.(6.3). Using these results in Eq.(6.19) we 
obtain, 

x\ W '\p(e,z)- V, e = f E aI"V,,). (6.22) 

Let us notice now that the actual length of the combination |VF _1 VF'| is 
equal to one since in one step of the random walk only one letter is added 
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or subtracted according to the definitions made above. Finally, let W = e 
and ,since \e\ = 0, we obtain from Eq.(6.22) the following result: 

E> 1 

A (6.23) 



Pz 
since 

=2n (6.24) 

by construction. Here P = P(e, z).By combining Eq.(6.13) with Eq.(6.23) we 
obtain 

A = r (6.25) 

n + (n 2 — (2n — l)z 2 ) 2 

Obtained result coincides exactly with that obtained by Kasteleyn,Ref[23],e.g.see 
his Eq.8 , where it was given without derivation. By combining Eqs(6.21) and 
(6.25) we obtain the final result for P(W, z) which can be used now for cal- 
culation of II(W).In particular, for W = e we obtain,using Eq.(6.20), 

n(e) = H(0) = -^— x (6.26) 

in complete accord with earlier obtained Eq.(6.16) . It should become clear 
by now,that consideration of random walks on groups provides practically 
as much information as one could possibly get by much more tedious calcu- 
lations using the real-space approach. This is so because the group-theoretic 
methods provide an equivalent and complete description of the underlying 
geometric and topological properties of a given manifold[3, 75-77] .We would 
like now to elaborate this observation in the next section. 



VII. Some applications 

The recurrence/transience of random walks has actually numerous physi- 
cal applications which had been mentioned already in our previous works [14]. 
To make our presentation self-contained, we would like to remind to our 
readers about some of these applications and to discuss additional physical 
applications in this section. 
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Figure 7: Planar projection of the trefoil knot 



A. Polymers 

First, the results of previous section can be used in the theory of poly- 
mer solutions [14] since they provide the exact mathematical formulation of 
the concept of entanglement through use of mathematically well defined 
concepts of recurrence and transience. It should be clear by now that the 
random walks on groups are not at all limited to the planar problems 
since the homotopy of paths are not limited to two dimensions. In partic- 
ular, one can think now about random walks in the presence of knots and 
links(the possibility mentioned already in our earlier work, Ref[14]). From 
this reference we know that almost all flexible polymers are either knotted 
or quasiknotted in solution in case if their length L is large. Formation of 
knots/links may affect both static and dynamic properties of polymer solu- 
tions and also polymer interfaces as is well known. Nevertheless, to detect 
the presence of knots in solution is very nontrivial problem. Fortunately, 
using methods developed in this work this problem may be solved to some 
extent. Let us explain briefly how this can be done on examples of trefoil 
,Fig.7, and figure eight ,Fig.8, knots. Using either results of our Appendix 
A2 of Ref[14] or more specialized books on knot theorye.g. Ref[43],[81], one 
obtains the following group presentations 

Gt =< a,b,c\cb = ba = ac> (7-1) 
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Figure 8: Planar projection of the figure 8 knot, 
for the trefoil and 

G 8 =< a,b\ > (7.2) 

for the figure 8 knot complement groups. We would like to remind our read- 
ers that the above group presentations are for the fundamental groups of 
the complements of these knots in S 3 (i.e.for S 3 \.K" where K is the knot 
of given type and S 3 = R 3 U {oo}).It is being assumed that the knot has 
finite thickness (framing) so that the toroidal-like tunnel occupied by K is 
removed from B? thus converting it to space with possibly very complicated 
topology .The Cayley graph for the fundamental group of the trefoil is de- 
picted in Fig.9 while that for the figure 8 is depicted in Fig.6. Fig.9 has 
rather complicated structure which was discovered by Max Dehn [82] whose 
Dehn's twists we also going to discuss below. When viewed from above ,one 
can see the three- valent Bethe lattice (the terminology used in physics liter- 
ature for the Cayley graph [74]) . At each vertex of such lattice three infinite 
ladders are glued together in an obvious way. To describe random walks 
on such lattices is not an easy task and, besides, the presentation given by 
Eq.(7.1) is not unique, for example,using results of Ref[81] one obtains as 
well. 

G T =<x,y\x 2 =y 3 = 1> (7.3) 

with the corresponding Cayley graph (part of it) depicted in Fig. 10. Clearly, 
the random walk in the complement of the figure 8 knot is transient (that is 
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Figure 10: An alternative universal covering space for the complement of 
trefoil knot. 
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Figure 11: Base space for the graph depicted in Fig. 10 



if one could make a solution of linear and circular knotted polymers which 
all are figure 8 knots ,the linear polymers will become entangled with such 
knotted structures so that it will become practically impossible to separate 
these two types of polymers). To determine the recurrence/transience for 
the complement of the trefoil knot requires some additional efforts which 
are worth spending,see the next subsection. Already now few remarks are in 
place. Just by looking at Fig. 10 one can make the following observations. If 
we collapse the solid line loops to lines and the dashed triangles to points, we 
would obtain three-valent Bethe lattice as depicted in Fig.9. This lattice is 
the universal covering space for the triangle depicted in Fig.ll. Actually, to 
be consistent with the existing conventions of combinatorial group theory [83] 
one has to modify Fig.ll as depicted in Fig. 12. This is needed because 
, according to the existing rules , at given vertex only two edges of the same 
color(e.g. solid lines and/or dashed lines) can meet. By looking at Fig.10 
this fact then explains why figure 8 complement is transient while that for 
the trefoil is recurrent. Indeed, for the figure 8 we have at each vertex the 
probability 3/4 to go forward and only 1/4 to go backward while for the 
trefoil we get 1/2 to go forward and 1/2 to go backward ,much like for the 
random walks on the line which represent the universal covering space for 
once punctured plane discussed in section 2. Surelysuch walks are recurrent 
as it had been demonstrated already. 

Since the entanglements affect both static and dynamic properties of 
polymer solutions,it is clear, that the presence of knotted structures should 
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Figure 12: The same base space as in Fig. 11 but with better resolution 
required by the graphical representation of groups [] 



be important in changing these properties. The presence of knotted struc- 
tures is also important biologially since for different knots the first return 
and the escape times should be different .If the diffusion processes take place 
in the vicinity of such knotted structures, e.g. DNA's, the outcomes will be 
different for different types of knots. This observation could be important in 
processes which involve the molecular recognition. Let us however return to 
the problem of transience/recurrence for the random walks on fundamental 
groups for thefoil and figure 8 in order to discuss additional very important 
aspects of the recurrence-transience problem. 

B. Classical and quantum chaos and random walks on Te- 
icmiiller modular group 

From the knot theory textbooks, e.g. see Ref[43,81], it is well known that 
although figure 8 and trefoil knots look very similar, they actually belong to 
completely different "universality classes" as it can be seen already from the 
difference in their group presentations and will be explained further below. 
The trefoil knot is just a typical representative of the class of torus knots 
while the figure 8 knot belongs to the class of hyperbolic knots (actually it 
is the simplest example of hyperbolic- type knots [84]). It is not our purpose 
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Figure 13: Topological structure of once punctured torus 



to go into all details of the above distinction between these different classes 
of knots. Instead, we are going to emphasize aspects of these knots which 
are logically related to the discussions we had presented so far. 

In our earlier work, Ref [1] , we had noticed the fact that the fundamental 
group of the the the punctured torus is the same as for figure 8. This could 
be easily seen from the following picture, Fig. 13. Although in physics litera- 
ture quantum mechanics of such "leaky torus" had been discussed [85] along 
with many generalizations [9- 12], including those which involve the quan- 
tum Hall effect [86] ,e.g.Aharonov-Bohm effect under the periodic boundary 
conditions, etc., surprisingly, the spectral analysis and other related aspects 
of the punctured torus are still at the forefront of research in mathemat- 
ics [87]. Because of this, we would like to emphasize here some aspects of the 
problems related to punctured torus which had not been (to our knowl- 
edge) discussed in physics literature. Following Ref. [88] ,it is convenient to 
group the generators a and b of the free group Gs according to the following 
scheme: 

ab — > L 

ba — > L 
(bb)- 1 -► L 
ab- 1 -► R 

bb -> R 
b- l a -> R 

where R and L stand for "right" and "left" as it is usually done in symbolic 
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dynamics [89]. Surely, the above correspondence (automorphism) is quite ar- 
bitrary and it is obvious that other possibilities could be tested. What is 
less trivial is to prove [90] that the totality of words generated by sequences 
given by Eq.(3.4) coincides exactly with the totality of words of the type 

w = R n 1L n 2 ( 74 ) 

where [n!,n 2 ,...] is continuous fraction expansion of some number ,say 9 < 1. 
Recall [91], that the continuous fraction for is actually given by 

9 = ^ (7.5) 

n\ H y 

n 2 + 



n 3 + ... 



The continuous fraction always stops if 9 is rational number or it becomes 
periodic for some irrational numbers, e.g. those which come from solutions of 
quadratic equations, and it is infinite and nonperiodic for generic irrational 
9 . We would like now to argue that for the torus-type knots 9 is always 
rational number while for hyperbolic knots ,e.g. for figure 8, 9 should be 
irrational . Since in the first case words have finite length this means that if 
we have one of such longest words ,then multiplication by R (or by L) should 
make such word trivial, i.e. W=l. The trivial word, when depicted graph- 
ically [83], simply means that there is a closed loop in the graph. This, in 
turn, means that the motion is recurrent .If there are no closed loops, then 
the motion could be either transient or recurrent. . We had established 
this fact already in the previous section for the figure 8 . Since the trefoil 
knot is associated with three- valent Bethe lattice it means that the random 
walks on such lattices behave quite differently as compared with Bethe lat- 
tices with vertices of even valency . This difference also should show up 
in various properties of statistical mechanical models on such graphs [80]. In 
this work , naturally, we are not going to touch these subjects. Let us now 
explain better how all this is associated with chaos. Superficiallyan infinite 
word, Eq.(7.4), made of totally random sequences of R and L already sug- 
gests chaoticity. In general, the word W may or may not be random . It 
may correspond to the sequence of base pairs along DNA, to some number 
written in binary system, the sequence of "up" and "down" spins in Ising 
-like model, to sequence of hits in the Poincare return maps, etc .One can 
do still better than just these observations. To this purpose some knowl- 
edge of the knot theory is helpful and we refer our readers to the specialized 
literature, e.g. see Refs[43,81]. 
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Let is begin with a very well known fact from the combinatorial group 
theory which is this. Let G be free nonabelian group made of two generators 
a and b. Let us construct a commutator aba~ 1 b~ 1 = /C, then the quotient 
G/JC is just an abelian group that is , instead of presentation given by 
Eq.(3.5), we obtain now 



and this is just the fundamental group of the torus. Such abelianization had 
effectively removed the puncture from the torus so that it became essentially 
the product of two circles with the group Z©Z of covering transformations 
instead of just Z considered in section 2. This procedure of abelianization is 
characteristic not only for the punctured torus but for Riemann surface of 
any genus g :it is sufficient to make just one hole in such surface in order 
to obtain the free group of 2g generators. Using this fact ,one can prove the 
following related theorem 

Theorem 7.1. // F is finite graph with ol\ vertices and 012 edges then, 
the fundamental group of the graph ~K\(F) is just a free group of n generators 
where n = 1 + a.2 — ot\- 

Proof.Please,consult Ref [43] ,page228. 

The above theorem and the results presented in previous sections ex- 
plain at once the results of Ref [13] where the hyperbolicity of the disordered 
network was established using much longer chain of arguments. 

Going back to the torus case, one might think that the quantum mechan- 
ics on the torus is trivial. To a large extent this is the case, as compared 
with the case of the punctured torus, and ,indeed, many results for toroidal 
topology had been obtained in the theory of dynamical systems[16], quantum 
billiards [92] and the conformal field theories [93]. Nevertheless, even this case 
is still under detailed investigation as recent literature indicates [94]. In any 
event, by analogy with the circle, we shall utilize the properties of the cov- 
ering space which is just the infinite two dimensional lattice. We would like 
to construct explicitly the group representation which describe the motions 
on such lattice(s). Fortunately, this has been done already so that we can 
use known results. Moreover,these results will be useful also for the punc- 
tured torus case. Since the torus T=S' 1 x S 1 and the case of S 1 was studied 
in section 2 ,we can begin our discussion by dissecting our torus so that 
it becomes a square(without loss of generality)with sides of unit length . 
Then, the location of some point on the torus(or inside the square) can be 
described in terms of the four real coordinates 



G a =< a, b \ab = ba > 



(7.6) 



(xi, X2, X3, X4) = (cos a, sin a, cos 0, sin (5) 



(7.7) 
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or, alternatively, in terms of two complex coordinates z and w such that 




(7.8) 



2 



2 = 1. But Eq.(7.8) is the equation for 3-sphere 
S 3 = {{z,w} € C 2 : \z\ 2 + \w\ 2 = 2} 



since \z 



= \w 



(7.9) 



It can be proven [95] (and the proof is not too difficult) that all torus knots are 
obtained as result of intersection of the sphere S 3 with the complex surface 
z v _|_ w i = q where p and q are relative primes so that p/q is some rational 
number. Although all this is mathematically correct, we would like to provide 
intuitively more convincing picture of what is actually taking place. For this 
reason, we need to use the universal covering surface which is square lattice 
with the side length being 2ir .For construction of the trefoil knot, we need 
to consider the basic rectangle with sides of length 47r and 6ir respectively 
which can cover the entire complex plane by obvious translations. Let us 
now draw the straight line from the point (0,0) to the point (4ir, 6ir) and let 
us by means of obvious translations bring the squares affected by this line 
to the initial square which corner is located at (0,0). All this is depicted on 
Fig. 14. Now, the pattern in this initial square can be transferred back from 
the universal cover to the dissected torus without any change. Actuallythis 
pattern could be created directly in the base space too as depicted in Fig. 15. 
This figure was created with account of the periodic boundary conditions on 
the torus and, already from this picture one can recognize the knot projection 
if one properly resolves the intersections between the dashed and the solid 
lines(as it is always required in knot theory [43]). Fig. 15 allows us to reglue 
the square into the torus knowing that , indeed, the three dimensional knot 
is going to be created as result of such operation. This is depicted in Fig. 16. 
On this figure two trefoils are depicted :one is (3,2) and another is (2,3) 
which are, actually, the mirror images of each other and are called the right 
and the left hand trefoils respectively [43]. Now, what all this has to do with 
our original intentions?. For this,we have to talk about the analytical theory 
of toral automorphisms.lt is clear,even without any actual calculations,that 
the figure 8 knot should be associated with the lines of irrational slope in 
the covering lattice, Fig. 14. Such plausible reasoning would be to a large 
extent correct if we would forget the results of section 3 in which we had 
established that the universal covering space of the punctured torus is not 
a square lattice but rather the Poincare upper half plane H 2 . This fact 
provides yet another evidence of profound distinction between the trefoil 
and figure 8 knots. The situation can be repaired somehow if we provide 
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Figure 15: Dynamics of the rational paths in the base space of the torus 
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different interpretation of the results depicted in Fig. 14. To this purpose 
one needs to introduce the notion of the Teichmuller space of the punctured 
torus [96]. To make this concept readily comprehensible let us start with 
the notion of Dehn twists. As it had been proven by Nielsen [74], following 
works of Dehn[82],all selfhomeomorphisms of two dimensional surfaces can 
be performed with help of sequence of Dehn twists. Fig. 17a depicts one of 
such twists, while Fig. 17b demonstrates what such twist can produce on the 
surface of a torus. This simple observation can now be broadly generalized. 
To do this several steps are required. To begin, let us consider ,instead of 
a square lattice, Fig. 14, more general lattice ,as depicted in Fig. 18, which is 
characterized by the complex parameter r(which is the ratio of periods oj\ 
and t02 (some relatively prime integers) along the respective axes) so that ,in 
accord with Eq.(7.6), the abelian group of translations F T can be described 
as 

r T = {7 = m + riT\m,n £ Z;t G C,Imr > 0} (7-10) 

Surely,the difference between different tori lies in different values of r. Let 
us reglue now the basic parallelogram in order to form torus with canonical 
marking,Fig.l8b. For the reference torus we can take,say r = i.Now,if we 
change r we thus going to form new torus which , evidently, cannot be 
isometrically transformed to the old one. The relation between the new 
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Figure 19: Fundamental domain for PSL(2,Z) 



torus and the old one is established through 

/ (it + b 

T = ; 

ct + a 



(7.11) 



with a, b, c and d being integers subject to constraint ad — bc = l.The group 
PSL(2,Z) of the transformations just described is evidently the subgroup of 
PSL(2,R) introduced in section 3 and is known as modular group. Since 
PSL(2,R) is the group of isometries of H 2 , evidently, PSL(2,Z) is also acting 
on H 2 . The moduli space M=H 2 /PSX(2, Z) for the torus, by definition 
, coincides with the fundamental domain , Fig.19, for the group PSL(2,Z) 
as is well known [51]. The Teichmuller space is related to M but is more 
physically appealing as it will become clear momentarily. To this purpose, let 
us discuss the sequence of transformations as depicted in Fig. 18c and d.From 
this figure it is clear, that upon regluing of the deformed parallelogram the 
pattern of canonical markings ,Fig.l8b, had undergone some changes and,in 
view of Fig. 17, these changes could be interpreted in terms of the sequence 
of Dehn twists made with respect to canonical markings. Let us discuss an 
example of the trefoil knot. Let we have originally the square lattice with 
sides of the basic square 1 and i respectively. To construct a new lattice we 
use as a guide Fig. 18. The fundamental paralleogram of new periods can be 
chosen according to the following set of equations written with respect to 
the "old" axes: 
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Figure 20: An example of an area preserving toral homeomorphism related 
to formation of the trefoil knot 



/(t) = loit = ai + b 



(7.12) 



where in Eq.(7.12) we have to put a = 3 and b = 2 in order to obtain the 
trefoil knot. The numbers 3 and 2 are precisely the numbers of full twists 
with respect to canonical basis. In addition to Eq.(7.12) we also have 



The above two additional equations reflect the fact that the transforma- 
tion leading from the original square to the final parallelogram are area- 
preserving which is insured by Eq.(7.14) while Eq.(7.13) comes as result 
of such constraint. With a = 3 and 6 = 2 Eqs. (7.14) produce the fol- 
lowing set of solutions for c and cZ:a)c=l,d=l or b)c=-2 and d=-l.Fig.s 20 
and 21 demonstrate all this graphically. Both solutions make sense for the 
following reasons. These solutions provide the number of Dehn twists for 
the second basic curve involved in canonical marking. This second basic 
curve intersects the first one in just one point ,Fig.l8b,initially and Fig. 



/(I) = coi = ci + d 



(7.13) 



and 



ad — be = 1 



(7.14) 
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Figure 21: Different toral homeomorphism associated with the same trefoil 
knot 



18d,finally.But the amount of twisting of this second curve is related to 
the first one only by the law of area conservation so that in both cases 
we get the trefoil knot for the first basic curve. One can express all this in 
terms of vectors (a,b) and (c,d) which provide the numbers of Dehn twists 
along "vertical" /"horizontal" directions for the first and second canonical 
curves respectively.After this discussion it becomes clear that if we reglue 
the parallelograms in Fig. 20 and 21, we reobtain respective toruses as de- 
picted in Fig.l8d. Evidently, if we physically relate the initial square to 
the final parallelogram(s) we recognize that some stretching (may be twist- 
ing also)was necessary in order to bring the square to the parallelogram 
before gluing. Hence, the Dehn twists are always associated with some 
stretching and twisting of the underlying surface. This fact is directly 
associated with the notion of the Teichmiiller space since it is just the 
space of parameters responsible for surface deformations. In our case, the 
driving parameter is r so that the Teicmiiller space is H 2 for both torus and 
punctured torus [96] as it should become now intuitively clear based on the 
discussion we had just completed. Finally, the sequence of Dehn twists asso- 
ciated with selfhomeomorphisms(diffeomorphisms) of torus (surface of genus 
g, in general) causes some fundamental changes in the marking responsible 
for the fact that, for example, different torus knots are isometrically differ- 
ent. Therefore, it is possible to associate the mapping class group, known 
also as Teichmiiller modular group, with the group generated by all 
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nontrivial (that is not reducible to identity) Dehn twists (associated with 
canonical marking) so that the random walk problem of section 6 becomes 
the problem about the random walk on the mapping class group (except 
,now we have to deal with random walks on groups which are not free). An 
excellent and readable description of the mapping class groups could be 
found in Ref.[97] while good discussion of random walks on various groups 
which are not free could be found in Ref[21].Our first illustration of nontriv- 
iality of this group is the isometric difference between the tori depicted in 
Figs. 20 and 21. Let us discuss this topic a bit further since it has some appli- 
cations to dynamics of liquid crystals and 2+1 quantum gravity as described 
in our earlier work, Ref [17], (and ^specially, Appendix of this reference where 
connections of these dynamical problems with toral automorphisms (in the 
simplest case)are discussed) . 

If we select r in the fundamental domain, Fig. 19, i.e. in the moduli space 
M, then, it can be shown[51] ,that different r £ M are not connected with 
each other by means of Mobius-type transformation, Eq.(7.11).But,once we 
picked some specific r ,then, this r is connected with other r' which lies 
outside the moduli space via transformation of the type given by Eq.(7.11) 
. By means of successive applications of transformation law Eq.,(7.11), 
one can reach all conjugate points in the corresponding hyperbolic triangles 
which tessellate H 2 so that application of Eq.(7.11) to all points of the moduli 
space will cover the entire H 2 . Thus the Teicmiiller modular group is acting 
on the moduli space. Moreover, it is possible to show[51] that the arbitrary 
matrix 



V = 



a b 
c d 



) 




(7.15) 



where 




(7.16) 



and 




(7.17) 



Since 




(7.18) 
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projectively,e.g.see Eq.(7.11), this explains why in expansion (7.15) all T's 
have the first power. The fundamental domain for PSL(2,Z) in Fig. 19 is 
depicted with generators T and S in mind. As for the generator S ,we had 
already encountered it ,e.g.see Eq.(3.6). This means :a) that the numbers 
n l5 ....,n/ can be arbitrary integers and b) that using Eqs(7.15)-(7.17) we can 
always replace a -generator in Eq.(3.6) with the appropriate sequence of T 
and S generators. Moreover,the combination TS = Q produces projectively 

Q 3 = 1 (7.19) 

so that,in view of Eq.s(7.3) ,(7.18) and (7.19), we can consider the random 
walk on the complement of the trefoil using the generators of PSL(2,Z). This 
result is in accord with Ref[81] where other arguments were used to reach 
the same conclusions. From the discussion we had just presented it is clear 
that the description of the toral automorphisms by means of the random 
walks on the mapping class group should be equivalent to traditionally used 
in the theory of dynamical systems [16]. E.g., for example, the torus mapping 
depicted in Fig.20 could be equivalently described in terms of the map given 
by 

y' = x + 2y (7.20a) 
x' = x + 3y (7.20b) 

and,of course, 

1-3-1-2 = 1 

The above arguments (x and y) in this map should be taken modi [16] by 
analogy with S 1 case earlier discussed in section 2. For the case of toral 
automorphisms the use of Teicmiiller space interpretation appears as plau- 
sible curiosity. The things change dramatically if the self homeomorphisms 
of surfaces of higher genus is of interest, e.g. in gravity or liquid crys- 
tals^]. In mathematics literature ,the Nielsen-Thurston theory of surface 
selfhomeomorphisms is well developed [24,98] and, in view of the discussion 
already presented, can be used for description of various physical phenom- 
ena, including dynamics of liquid crystals and gravity .Indeed,based on 
the results just obtained the toral selfhomeomorphisms could be of three 
types:a)reducible,when the words like that in Eq.(7.15) are of finite length 
(this corresponds to finite continued fraction expansion,Eq.(7.5));b)periodic(which 
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may correspond to the periodic continued fraction of some quadratic ir- 
rationals)and ,c) irreducible (or Anosov (pseudo Anosov in more general 
case)) which may be associated with truly irrational continued fraction ex- 
pansions. The first possibility in the language of liquid crystals[99] corre- 
sponds to solid phase,the second to the hexatic phase and the third to the 
liquid phase. For more details,please, consult our earlier work,Ref[17]. 

C. Three-manifolds which fiber over the circle 

Connections between theory of dynamical systems and theory of knots 
and links are the most naturally described in terms of the Nielsen- Thurston 
theory[98] and , although this field is still under active development even 
in mathematics [100] , we would like to provide here a sketch of what one 
may expect upon development of such connections. This information may 
be useful in fields ranging from gravity through liquid crystals, dynamics of 
fracture and even folding of proteins (for references related to these problems 
, please, consult our earlier work,Ref[17]). 

Description of dynamical systems in terms of Lyapunov exponents, frac- 
tal and multifractal dimensions,entropy, etc is well known[101]. Recently, 
however, more and more emphasis is being made on connections between dy- 
namical systems and the theory of knots and links [102] . The complements 
of knots and links are 3-manifolds and, moreover, there are 3- manifolds 
which are not related to links /knots ^specially those which are not com- 
pact. The 3-manifolds arise completely naturally in the theory of dynamical 
systems. To realize this fact it is sufficient just to take a look again at Fig. 15. 
Extension this picture to a cube with periodic boundary condition creates 
already our first 3 manifold, i.e. cube with the opposite sides properly iden- 
tified. This example already provides a clue as to how 3-manifolds can be 
constructed: e.g. take some polyhedron and identify the opposite sides in 
some systematic way. As a result of such an identification we again obtain 
some 3-manifold. Since some of these 3 manifolds are related to knots it 
is interesting to know that, for example, Fig.8 knot complement ,which is 
also a 3-manifold, can be obtained by proper gluing of faces of two ideal 
tetrahedra to each other. This and other very beautiful illustrations of 
making 3-manifolds associated with simple knots/links could be found in 
Ref.[103]. In case of trefoil ,the associated 3-manifold is a (3.2) Lens space. 
It is easily created out of two solid toruses.To do so we should have one 
torus with standard canonical marking being glued to other (twisted) torus 
which contains the trefoil knot on its surface so that the initial (canonical) 
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marking is glued to the final marking and then the rest of two torus sur- 
faces being glued to each other [104]. This example gives an idea about 
the way 3-manifolds are related to surface dynamics. Let us now be more 
specific. Let S be some surface and let xG S .Consider a homeomorphism 
ijj : S — > S . Imagine now that this homeomorphism takes place in stages 
,i.e.evolvs in time so that at time t=0 ip(x) = x while at time 1 tp(x) may 
or may not coincide with x Vx G S. The time 1 is, of course, quite arbitrary 
and is chosen for convenience only.Let us now identify (x,0) with (ip(x), 1) 
Vx G S'.This identification leads us to a 3-manifold T^, . This 3-manifold 
is a fiber bundle with fiber S, base S 1 (a circle) and monodromy ip. It is 
clear that the monodromy ip is directly related to the random walk on the 
Teichmiiller modular group (as described in the previous subsection). It is 
clear also that different surface dynamics will lead to different 3-manifolds 
in general so that the topology and geometry of 3-manifolds is yet another 
way of description of dynamical systems [100, 102]. 

The above picture can be elaborated now as follows. Since there is one-to- 
one correspondence between the sequence of Dehn twists and the sequence 
of words being regulated by the continued fractions expansions it is clear 
that the trivial fiber bundle can be associated with the reducible selfhome- 
omorphisms and that the periodic selfhomeomorphisms produce naturally 
3-manifolds known also as Seifert fibered spaces [8 1,84, 104]. Surely, all torus 
knots are associated with periodic selfhomeomorphims and,therefore are in 
one- to-one correspondence with Sefert fibered spaces as was already estab- 
lished by Sefert in 1933(the sufficiently simple proof of this could be found 
in chr.6 of Ref.[81]).What is much less clear why the above construction of 
3-manifolds should be also applicable to irreducible (or pseudo Anosov) 
selfhomeomorphisms. The fact that this is actually possible had been proven 
rather recently [100]. And ,the most remarkable fact of this proof lies in use 
of Feigenbaum-Cvitanovich results for one dimensional maps which, in the 
light of the results we had presented already should not be totally unex- 
pected. 3-manifolds associated with knots such as figure 8 are hyperbolic. So 
that all knots other than torus knots should be associated with hyperbolic 3- 
manifolds. That this is indeed the case had been proven by Thurston[86,103] 
and was mentioned already in our earlier work, Ref [14] . The hyperbolicity of 
knots is directly connected with transience of random walks in the presence 
of such knots as we had already discussed in section 7. A. Moreover, using 
the results of hyperbolic geometry it is possible to get information about 
fractal dimensions of the limits sets (in case of dynamical interpretation). 
Some introduction to this subject can be found in our earlier work,Ref.[105]. 
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D. Schottky doubles, Markov triples and Grothendieck's 
dessins d'enfants 

Although we had covered many aspects of the figure 8 story, there are 
many topics which we had not mentioned at all so far. In this subsection 
we would like to correct this deficiency. Evidentlyour presentation is going 
to be brief but, since the topics to be discussed are directly and logically 
related to the rest of this paper, it is impossible not to mention about them 

Let us begin with the notion of Schottky double. From the discussion we 
had so far it is clear that treatment of boundaries in the hyperbolic geometry 
is not as simple as it is in Euclidean case. For instance, suppose we would like 
to calculate the spectrum of Laplacian on the punctured torus .The question 
immediately arises: how the spectrum of Laplacian will depend upon the 
boundary conditions at the puncture ? To some extent this problem was 
solved by Gutzwiller[85] who however admits, e.g. see page 376 of Ref.[85], 
that more systematic mathematical treatment is given in the monograph 
by Lax and Phillips [106]. Surprisingly, the case is still not closed as recent 
Annals of Mathematics paper [87] indicates. Hence, it would be too naive 
to address this problem in its full generality in this paper. Nevertheless, we 
would like to mention one simple case which has its familiar analogue, for 
example, in standard quantum mechanics. The idea belongs to Schottky 
and is described in Ref. [96]. Consider a Riemann surface R with punctures 
and its copy R obtained by reflecting R in the mirror. Given these two 
copies , glue them along the boundaries of punctures . Thus obtained new 
Riemann surface is known as Schottky double. In the case of punctured 
torus the Schottky double is a double torus which ,not surprisingly, was 
also considered by Gutzwiller[85]. Such construction is equivalent of having 
perfectly reflecting boundary conditions in standard quantum mechanics and 
also has some applications to polymer solutions in confined geometries [14] 
, etc. The boundary conditions (that is the existence of puncture), as pointed 
by Gutzwiller [85] and known from the theory of Riemann surfaces[107] 
, could be modelled by imposing certain constraints on the commutator K, 
defined before Eq.(7.6).In particular, one can require 

K =aba- 1 b- 1 = V (7.21) 

(where V is transformation of the type given by Eq.(7.11)) to be parabolic 
, that is tr(P)=2. Such requirement provides one constraint on the explicit 
form of matrices a and b. Additional constraints could be imposed by locating 
vertices of the hyperbolic triangle at points 0, 1 and oo as depicted in Fig.2, 
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or hyperbolic quadrangle at points -1,0,1 and +00 , as depicted in Fig.3. 
Let us check this result explicitly. Using results of section 3 we obtain, 



a(z) 



z 



(7.22) 



2z+ 1 



b(z) 



z + 2 



(7.23) 



Accordingly, 




z 



(7.24) 



-2z + l 



and 



b(z)- 1 =z-2 



(7.25) 



Using Eqs.(7.22)-(7.25) in (7.21) a short computation produces 



-z 4- 1 

2^ "1" 8 



(7.26) 



From here,tr("P(z))=2 as required. The above simple calculation teaches 
us some very important lessons. First , since we are dealing with projective 
transformations ,e.g. that given by Eq.(7.11), we can always normalize our 
matrices the way which is the most convenient for us. In particular, al- 
ready in more than hundred years ago Fricke[108] had shown that for any 
unimodular matrices A and B 



(trA) 2 + (trB) 2 + (trAB) 2 = (trA) (trB) (trAB) + trK+2 (7.27) 



where, as before, K, is the commutator of A and B. Now , let tr/C = —2. This 
is quite permissible in view of Eq.(7.26) since multiplication by -1 of both 
numerator and denominator changes nothing(e.g.see Eq. (7. 18)). Such nor- 
malization however, makes a difference in Eq.(7.27) which now could be 
brought to form 



where C= (AB) _1 and the fact that tr(D)=tr(D _1 ) for any unimodular ma- 
trix D was used. By setting \trA\=3mi, \trB\=3m,2 and \trC\=3rri3 we obtain 
, instead of Eq.(7.28),the following equation for the Markov triples[25] 



(trA) 2 + (trB) 2 + (trC) 2 = (trA) (trB) (trC) 



(7.28) 



mi + TTI2 + m3 = 317117712^3 



(7.29) 
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with mi,m2,m,3 being nonnegative integers. By looking at this equation we 
immediately recognize that it is trivially satisfied by mi = 777-2 = ^3 = 1- 
This information ,when combined with the requirement of unimodularity, is 
sufficient for restoring the explicit form of matrices A and B ,and , hence, 
AB. Questions arise :a)are there other solutions to the equation for Markov 
triples and b) what is their physical significance if they indeed exist? The 
answer to the first question is known in mathematical literature and will 
be discussed shortly below. The answer to the second question could be 
related to the fractional quantum Hall effect. Indeed, if solutions other than 
trivial to Eq.(7.29) do exist, ,then, based on the results of sections 3 and 
4, we conclude that they may be associated with different shapes of the 
fundamental domains in H 2 plane and ,accordingly,to different locations of 
punctures in S 2 .Thus if the punctures are associated with anyons, then , 
evidently, their mutual position cannot be arbitrary so that we have some 
sort of solid phase in the cluster approximation (made of 3 punctures). 
Within this approximation the filling fraction[14] ,0<^ < 1 can be calculated 
and, because the pattern or Markov triples is very intricate (see discussion 
below) one obtains a very rich array of the permissible filling fractions char- 
acteristic for the fractional quantum Hall effect[109].To what extent such 
cluster picture can be extended to more than 3 punctures is related to the 
works of Grotendieck on dessin d'enfants[22] to be also discussed below. 

Looking at Eq.(7.29) it is clear that there is a complete symmetry be- 
tween mi, 7772 and 7773. This symmetry can be broken if we notice that (2,1,1) 
is also solution to Eq.(7.29). This solution is called singular [25]. It could be 
shown ,that this is the only solution for which 2 out of three Markov numbers 
are the same. All other solutions have distinct Markov numbers. To generate 
these other solutions (Markov triples) it is sufficient to have a "seed" in the 
form (1771,7772,7773) which ,by definition, obeys Markov equation, Eq.(7.29). 
Then, the first generation of new Markov triples given by 

(7771,7772,7773), (mi,?77 2 ,?773) and (7771,7772,7773) 

where 

/ / / 

777i = 377727773 — 777i ; 777 2 = 3777l?773 — 7772 5 777 3 = 3?77l?772 — 7773 (7.30) 

.If one starts with (1,1,1) then, the Markov tree (its part of course) is 
depicted in Fig.22. The emergence of numbers 777 1 ,etc. , can be easily 
understood if we consider , instead of Eq.(7.29), its equivalent form given by 

f(x) = X 2 — 377727773X + 7772 + m 3- (7-31) 
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Figure 22: Fragment of the Markov tree closest to the original "seed' 
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Evidently, f{x) = for x = mi and x = 3m2Wi3 — mi. Additional infor- 
mation about Markov triples could be found in Refs. [110-112]. Now the 
question arises as to what extent one can generalize these results to more 
than tree punctures. There are several ways to proceed. For example, prior 
to Theorem 7.1. we've noticed that it is sufficient to make one puncture 
(hole) in surface of any genus g in order to obtain a free group with 2g gen- 
erators. By analogy with Eq.(7.29), one may think of solutions of equation 

x\ + x\ + • • • + x 2 n = nx\X2 ■ ■ ■ x n (7.32) 

as it was done by Hurwitz[113] and elaborated by Hirzebruch and Za- 
gier[114]. There is , however, much more striking approach to the whole prob- 
lem . It was proposed by Grotendieck [22] and was elaborated by Belyi[115] 
and others [22]. These mathematical results had found already their place in 
physics literature in connection with problems related to discretized string 
theory[22] and the theory of polygonal billiards [116]. As it was argued above, 
these results could be potentially useful in the theory of fractional Hall effect 
and/or quantum chaos, etc. 

The idea of the Grotendieck approach is rather simple. Take for ex- 
ample the torus case. The torus ,Fig.l8a(or 18b) can be made out of two 
triangles and accordingly, any Riemann surface of genus g also can be tri- 
angulated. Such triangulation of arbitrary closed Riemann surface could be 
easily understood if we would consider the motion of ficticious particle in a 
triangular billiard whose angles are some rational fractions of 7r. Instead of 
watching the complicated trajectory inside the triangle it is much better to 
unfold the trajectory in a way shown in Fig. 23. Such unfolded trajectory 
is just a straight line just as in the case of a torus ,Fig. 14. From the 
covering surface ,e.g. that depicted in Fig. 14, one needs now to restore the 
corresponding (Riemann) surface. This was done initially by Richens and 
Berry [11 7] and was subsequently generalized and elaborated by Zemlyakov 
and Katok[118] and many others,e.g.seeRef.[119]. Strictly speaking , to tri- 
angulate Riemann surface one needs to consider such triangle groups whose 
action on the fundamental domain for such groups will cover (tessellate) 
given Riemann surface without gaps or overlaps. Grotendieck calls such 
groups "cartographic groups". A good survey of the cartographic groups is 
given in Ref.[22].The triangulation of Riemann surfaces is not a new idea. The 
new idea lies in utilizing the results for the trice punctured sphere, sections 3- 
4, in order to achieve the triangulation of Riemann surface of any genus. The 
dessins d'enfants (the drawings of children) problem is directly related to 
the triangulation problem and can be formulated as problem of drawing an 
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Figure 23: Unfolding of the particle trajectory in triangular billiards: each 
time the particle hits the wall, one should extend the particle motion by 
considering its motion in the mirror image of the triangle with respect to 
this wall. 
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a) b) 

Figure 24: Triangulation of the arbitrary Riemann surface a) and the as- 
sociated with it double triangle b) which (upon regluing) is topologically 
equivalent to the trice punctured sphere as discussed in the text 



arbitrary(but pre assigned) closed graph on Riemann surface. This happens 
to be not a simple task .To make a connection between the trice punctured 
sphere and some Riemann surface one needs to have a mapping function 
which connects the triangulation on the Riemann surface with triangles on 
the sphere. This function had been constructed explicitly by Belyi[115] . We 
follow , however, the work of Shabat and Voevodsky[120] in order to illustrate 
the basic idea of the method. 

If the surface is orientable, it is possible to color the adjacent triangles 
in two different colors,say, black and white. Given such coloring, one can 
distinguish three types of vertices :o ,• and *.Let the red arrow connect o 
and*, let also the black arrow connect * with • and, finally, let the blue 
arrow connect • with o .Let us next associate the vertex o with number 0, * 
with 1 and • with oo . The above numbers are associated with the location 
of punctures on S 2 for the trice punctured sphere as discussed in sections 3 
and 4. To understand the connection, observation of Fig. 24a is helpful at this 
time. Let us select two adjacent triangles which ,by construction, are colored 
differently, e.g. see Fig.24b. Finally, if we glue together these triangles along 
the free edges, thus obtained object is just a trice punctured sphere. This can 
be easily understood if we imagine that the triangles are made of rubber and 
we blow up air through one of the open holes and then close it. Hence, to 
get a triangulation of an arbitrary Riemann surface 1Z , one shoold only 
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look for the appropriate branched covering with branches over the points 
0, 1 and oo. Clearly, for each point, say, on S 2 , there is a finite number 
of the corresponding points on 1Z , etc. The mapping function F(z) which 
actually maps S 2 \{0, 1, oo} to 1Z is given by 

( m + n ) m + n 

F{z) = [m+ y n z m {\ - zT (7.33) 
m m n n 

with m, n being some integers which value depends upon the genus of 1Z . 

From this discussion it is clear that thus obtained Markov triples could 
be now used for an arbitrary 1Z and this fact then proves their relevance for 
computation of the multitude of the filling fractions v characteristic for the 
factional Hall effect[109]. To complete our discussion , we have to consider 
still one more topic. 



E. Connections with braid groups 

We include this subsection in this paper in order to emphasize some 
aspects of theory of braids which are directly related to the discussion we 
had so far. We assume ,that the reader is familiar with some basic concepts 
of the theory of braids [81]. The braid group is infinite group of n-1 generators 
0"i,cj2, ...o n -\ obeying the following set of defining equations(in the planar 
case only) 

UiUj = OjOi for |i-j| > 2 (7.34) 

and 



(Ti(Ti +1 <Ji = a i+1 cTi(Ti +1 for 1 < i < n - 2. (7.35) 

As it was shown already by Artin (inventor of the Braid group) [121] , the 
braid group can actually be generated just by its two elements 

a = cria 2 ...(Jn-i (7.36) 

and 



<7 = ffi. (7.37) 

Thus, 

a i = a i - 1 aa- { - i - 1 \ (7.38) 
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and 



a n = (aa)" 1 - 1 . (7.39) 

In particular, let n=3,then, instead of Eqs(7.35), we get 

0"1<72<71 = <T 2 (TiCr 2 (7.40) 

Accordingly, 

a = o\G2 (7-41) 
and,using Eq.(7.40), it is convenient to introduce yet another generator 

b = aa (7.42) 



Then,Eq.(7.39) is reduced to 



a 3 = b 2 (7.43) 



Thus, we had obtained back the presentation group for the trefoil given ear- 
lier by Eq.(7.3). a n , introduced in Eq.(7.39), has some physical meaning . 
If we would have a trivial braid made of n strands , then application of a n 
to such trivial braid amounts of full rotation of the lover side of the frame 
while keeping the upper side fixed [81]. Therefore , if the system ,say, is not 
rotationally invariant, one cannot use the concept of braids while the con- 
cept of fundamental group of free generators still survives. The connections 
between braids and free groups could be made more explicit if needed . Here 
we would like only to touch upon the most immediate connections. To this 
purpose we need to introduce the notion of colored braid. Basically, the 
braid is colored if the end of each strand lies exactly over its beginning (that 
is for the ordered sequence on the bottom frame there is exactly the same 
ordered sequence on the top). From this definition, we notice that colored 
braids correspond to the distinguishable particles in physics terminology 
. For the colored braids P n made of n strands connection with free groups 
becomes very direct. Fig. 25 provides an example of a generic action of the 
braiding generator (an element of a free group) while Fig. 26 illustrates 
how the colored braid P4 could be made out of composition of braids made 
with help of the generators of the free groups F\ , F2 and F3 . For general n 
mathematically,such braid can be written as semidirect product 

P n = F 1 «(F 2 <(....(F n _ 1 «F n )) (7.44) 

For the colorless case, one can analyze the connections further by using 
Eq.s(7.36)-(7.39). 
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Pj p, 




(i) 

Figure 25: Action of the free group braiding generator a,- on the braid 




Figure 26: Colored braid made with help of composition of actions of gen- 
erators of free groups Fi,F2 and F3 made of one , two and three generators 
respectively 
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